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PREFACE 


■"PHE  present  work  is  intended  for  the  use  of  students  who 
are  learning  the  practice  of  Interpolation  and  Numerical 
Integration.  The  advantages  of  a  practical  kno wedge  of  this 
part  of  Mathematics  are  so  obvious  that  it  is  needless  to  insist 
on  them  here  :  and  these  subjects  form  an  important  part  of  the 
course  in  the  modern  Mathematical  Laboratory.  There  are, 
however,  so  many  claims  on  the  time  of  students  that  the  extent 
of  this  course,  as  of  all  others,  must  be  kept  within  narrow 
limits  :  and  it  has  therefore  been  necessary  to  restrict  the  treat- 
ment to  the  most  central  and  indispensable  theorems.  A  large 
number  of  numerical  illustrations  and  examples  has  been  given 
of  a  kind  likely  to  occur  in  the  applications  of  Mathematics. 

My  thanks  are  due  to  Professor  Whittaker  and  to  my 
colleague,  Mr  E.  M.  Horsburgh,  M.A.,  B.Sc,  Assoc.  M.Inst. C.E., 
for  their  valuable  criticisms  and  suggestions. 

D.  G. 

The  Mathematical  Laboratory, 

University  of  Edinburgh, 

July  1915. 
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PRELIMINARY    NOTE    ON    COMPUTATION. 

DEFORE    inti'0(hiciiig    the    formulae    required    in    interpolntion 
and  numerical  integration,   it  may  be  advisable   to   mention 
a  few  points  which  will  facilitate  the  work  of  computing. 

Where  computation  is  performed  to  any  considerable  extent, 
computer's  desks  will  be  found  useful.  Those  used  in  the  mathe- 
matical laboratory  of  the  University  of  Edinburgh  are  3'  0"  wide, 
r  9"  from  front  to  back,  and  2'  6i"  high.  They  contain  a  locker, 
in  which  computing  paper  can  be  kept  without  being  folded,  and 
a  cupboard  for  books,  and  are  fitted  with  a  strong  adjustable 
book-rest.  Thus  the  computer  can  command  a  large  space  and 
utilize  it  for  books,  papers,  drawing-board,  arithmometer,  or 
instruments.  Each  desk  is  supplied  with  a  copy  of  Barlow's  tables 
(which  give  the  square,  square  root,  cube,  cube  root,  and  the 
reciprocal  of  all  numbers  up  to  10,000),  a  copy  of  Crelle's  multi- 
plication table  (which  gives  at  sight  the  product  of  any  two 
numbers  each  less  than  1000),  and  with  tables  giving  the  values 
of  the  trigonometric  functions  and  logarithms.  These  may, 
of  course,  be  supplemented  by  a  slide  rule,  or  any  of  the  various 
calculating  machines  *  now  in  use,  and  such  books  of  tables  as 
bear  particularly  on  the  subject  in  question. 

Success  in  computation  depends  partly  on  the  proper  choice  of 
a  formula  and  partly  on  a  neat  and  methodical  arrangement  of  the 
work.  For  the  latter  computing  paper  is  essential.  A  convenient 
size  for  such  a  paper  is  26"  by  16";  this  should  be  divided 
by    faint  ruling    into    |"    squares,    each    of    which    is   capable    of 

*  For  descriptions  of  these  see  Modem  Instruments  and  Methods  of  Cal- 
culation, edited  by  E.  M.  Horsburgh.     London  :  G.  Bell  <fe  Sons.     1914, 
e.  1 
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holding  two  digits.  It  will  be  found  conducive  towards  accuracy 
and  speed  if,  instead  of  taking  down  a  number  one  digit  at  a  time, 
the  computer  takes  it  down  two  digits  at  a  time,  e.g.,  instead  of 
taking  down  the  individual  digits  2,  0,  4,  7,  6,  3,  it  will  be  found 
better  to  group  them  together,  as  20,  47,  63.  Every  computation 
should  be  performed  with  ink  in  preference  to  pencil;  this  not 
only  ensures  a  much  more  lasting  record  of  the  work  but  also 
prevents  eye-strain  and  fatigue. 

It  need  hardly  be  emphasized  that  seven-place  accuracy  cannot 
be  obtained  by  the  use  of  four  figure  tables.  A  fairly  safe  rule  is 
to  make  use  of  tables  containing  one  digit  more  than  the  accuracy 
required ;  to  use  more  accurate  tables  than  these  is  simply  to 
increase  the  amount  of  labour  without  increasing  the  elhciency. 
For  the  same  reason  contracted  methods  of  multiplication  and 
division  should  be  adopted  when  machines  and  tables  are  not  at 
hand. 


CHAPTER     I. 

SOME    THEOREMS     IN    THE    CALCULUS    OF 
FINITE     DIFFERENCES. 

1.  Differences. 

Let  -^/(x)  denote  the  increment  of  a  function  /(x)  correspond- 
ig  to  a  given  constant  increment  w  of  the  variable  x,  so  that 

A/{x)=/{x  +  w)-/(x). 

The  expression  A/{x),  which  is  usually  called  the  First  Differ- 
nce  of  the  function  f{x)  corresponding  to  the  constant  increment 
:;  of  the  variable  x,  will,  in  general,  be  another  function  of  x,  and 
s  such  will  have  a  first  difference  which  will  be  obtained  by 
perating  on  ^/(x)  with  A.  Calling  the  result  of  this  operation 
i-/{x),  we  have 

^\fix)  =  A{^/{x)} 

=  ^{/{x  +  w)-/(x)} 

=  {/{x+2tv)  -/{x  +  w)}  -  {/{x  +  zr)  -fix)} 

=/{x  +  2w)  -  2f{x  +  w)  +f{x). 

^his  is  termed  the  Second  Difference  of  f{x)  corresponding  to  the 
onstant  increment  w  of  x.  By  repeating  the  process  we  obtain 
n  turn  A-'y(x),  A^/{x),  ...  A" /(a;),  which  are  called  the  3rd, 
th,  ...  wth  differences  of  /{x)  corresponding  to  this  constant 
ncrement  w. 

Note. — It  is  always  possible  by  means  of  a  suitable  trans- 
ormation  to  make  the  value  of  w  unity. 

For  since  w  is  the  increment  of  x  in  f{x),  then  unity  is  the 

ncrement  of  —  .     Replacing  x  in  /{x)  by  wy,  we  obtain  a  new 

unction  F  {y),  which  is  such  that 

^F{y)  =  F{^J+\)-F{y). 
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Example  7.— Tabulate  and  difiference  the  successive  values   of  f{x)  =  x- 
from  a;=llO  to  x=118. 


X 

ar* 

A 

A2 

110 

1331000 

36631 

111 

1367631 

37297 

666 

112 

1404928 

37969 

672 

113 

1442897 

38647 

678 

114 

1481544 

39331 

684 

115 

1520875 

40021 

690 

116 

1560896 

40717 

696 

117 

1601613 

41419 

702 

118 

1643032 

A3 


A* 


In  the  above  table  the  column  headed  x^  was  obtained  from  Barlow's 
Tables.  The  column  A  gives  the  first  differences,  being  obtained  by  sub- 
tracting each  eatry  of  the  column  x^  from  the  entry  following  it,  e.g.  if 
a:=114,  a;+l  =  115,  A/(2:)  =  (a;  +  1)3- a;3=  1520875  -  1481544  =  39331.  This 
result  is  set  on  a  line  midway  between  the  numbers  corresponding  to  a;=114 
and  x=  115.  In  the  same  way  the  2nd  ditferences,  A-,  are  obtained  from  the 
1st  differences,  the  results  being  set  on  a  line  midway  between  the  latter, 
and  so  on. 

It  will  be  noticed  that  the  3rd  differences  of  x^  are  constant  and  that  its 
4th  differences  are  zero.  It  will  be  shown  later  that  the  7ith  differences  of  a 
function  of  the  nth.  degree  are  constant  and  that  the  differences  of  higher 
order  are  all  zero. 


Example  2. — To  determine  the  successive  differences  of  sin  {ax  +  b). 
Let    f{x)  =  si\i{ax  +  h). 
Then    A/(cc)  =  sin  {a  x  +  w  +  b)  -  s,\n{ax  +  h) 
=  2  sin  ^ a  w  cos  [ax  +  b  +  ^aiv) 


=  2 sin  \aiosin{ax  +  b  +  haw  +  ir) 
A'^fix)  =  A{  2  sin  ^  a  w  sin  {ax  +  b  +  i  a  w  +  tt)} 


-2sin  ^aw {sin  {ax  +  w  +  b  + -^aiv +  ir)-sin  {a x  +  b  +  l  aw +  Tr)} 
=  {2  sin  ^  a  iv}^  cos  {a  X  +  b  +  ^YowTtt) 


=  {2sin^aM;}2sin  {ax  +  b  +  ^2aw  +  2ir) 
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Similarly,  it  can  be  shown  that 


A^/(a;)  =  {2  8in  h  a  iu\^  sin  (a  x  +  b  +  ^  3  a  w  +  '^  n) 
ind,  more  generally,  that 

A"/(a;}  =  {  2  sin  ^  aw}"  sin  (ax  +  b  +  ^n  aio  +  7nr). 
This  can  be  proved  by  induction  in  the  following  way  :- 


l"+^/(x)  =  A  [{ 2  sin  i  a  w}"  ain  {ax  +b  +  ^n  aw +  mr)] 


{2sin^aio}n  [sin  {ax +  w  +  b  +  ^n  aw +  mr) 
-sin  {ax  +  b  +  ^naw  +  nir)] 


I  2  sin  ^  a  w  }''+i  cos  {ax  +  b  +  ii{n  +  l)aw  +  mr) 


=  {2sinia  ;/;}«+isin  {ax  +  b  +  i{n+ l)aw  +  {n+l)ir). 
But  the  theorem  has  been  proved  for  n—l,  2.     Hence  it  is  true  always. 
Example  3. — Find  the  first  differences  of  the  functions 

cos  {ax  +  b);  tan— i  ax  ;  a**. 
Example  4. — Find  the  ?ith  differences  of  the  functions 

—  ;    a*;   cos''(a.f +  6). 

x 

2.    Properties  of  the  Operator  A. 

The   operator    A    obeys   certain    of    the    fundamental    laws    of 
alsebra,  viz. — 

(i)  The  Law  of  Distribution. 

For  A  {/{x)  +  <^  {x) }  =  {/{x  +  iv)  +  4>  (x  +  w)}  -  {/{x)  +  i^  (x)} 

=  {  f{X  +  W)  -f{x))  +  {4>{x  +  w)-  c/)  (x) } 

=  A/(x)  +  A</,(x). 
(ii)  The  Law  of  ConwitUation  if  the  coefficients  of  the  various 
terms  are  constants. 

For,  if  a  is  a  constant, 

A  {a  .f{x)'f  =  af{x  +  w)  -  af{x) 
=  a{f{x  +  w)-f{x)} 
=  a.^f{x). 
(iii)  The  Laiu  of  Indices. 
For  A'- .  yf{x) 

=.  {A .  A  .  A  ...  (?■  factors;  .A    A  .  A  ...  {s  factors)} /(a;) 
=  {A  .  A  .  A  ..;  {r  +  s  factors)} /(a;) 
=  A'+'/(^). 
The    indices    r   and    s   are,  of    course,   integers:    the    operator 
A"  was  defined  in  the  previous  section  only  for  integral  values  of  ti. 
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Example  l.—lif{x)  is  a  polynomial  in  x  of  the  nth  degree,  then  its  nth 
difference  is  constant  and  its  differences  of  higher  order  than  n  are  zero. 

Let  /(x)  =  a„a,-" +«»-iX"-i+  ...  +ai  x  +  wo 

in  which  the  coefficients  a„,  a„_i,   . .  aj,  ao  are  constants. 

Then     ^f(x)  =  {a„  (x  +  /'■)"  +  «"-'  {x  +  w)"-'^  +  ■■■  +a,{x  +  v;)  +  a^  } 
—  [a,,  X"  +  a„_i  X"-'  +  ...  +  rti  a;  +  ao  } 
—  n  a„  a;"-'  c  +  /'„_o  x"  --  +  ...  +  ?*i  x  +  to 

where  6„_2,  ...,  '>i,  60  are  constants. 

We  thus  see  that  A/(a;)  is  a  polynomial  in  x  of  degree  (n  - 1).  Repeating 
the  operation,  we  have 

:^-f(x)  =  {n  a„  (x  +  ;y)"~^ '"  +  &n-2  (a^  +  "')"--  +  •  •  •   +K{x^  "')  +  ?'o} 
-{?ia„a;"-i  «;  +  &„_2a;"-2+  ...  +?)ja;  +  &o} 
=  n(w- I)  a„a;"-2  H;2  +  c„_3a;''-2+  ...  +Cia;  +  Co 
where  c,i_3,  ...,  Cj,  Co  are  constants. 

Hence  A-/(^)  is  a  polynomial  in  x  of  degree  (?;  -  2).  It  is  evident  that 
by  continued  repetition  of  this  process  the  degree  of  the  function  will  be 
reduced  to  zero,  and  that  we  shall  obtain 

Lnf(x)  =  ?t  (h  -  1 )  (vi  -  2)  . . .  3  .  -2  .  1  a„  )r«  , 

i.e.,  the  nth  difference  of  the  function  is  a  constant.  Obviously  the  next  and 
all  differences  of  higher  order  will  be  zero.  This  result  is  of  greo,t  import- 
ance ;  probably  no  other  theorem  will  be  iipplied  as  frequently  a.s  this.  For 
instance,  it  enables  us  to  complete  a  mathematical  table  in  which  there  are 
gaps.     This  is  illustrated  in  the  next  example. 

Example  2. — If 
/(x)=:0-5563,  /(cef  I)  =  0-5682, /(a:  +  3)  =  0-5911,  /(a;  +  5)  =  0-6128, 

obtain /(a: +  2)  and /(a; +  4)  on  the  assumption  that /(a;)  is  a  polynomial  of 
the  third  degree. 

It  has  been  shown  in  §  1  that 

AV(a--)  =/(■^•  +  2  w)  -  2f{x  +  w)  +f(x). 
Bj'  repeating  the  process  we  get 

AV (a;)  =f(x  +  A  i.)  -  4/ (a;  +  3  n-)  +  6/(x  +  2  *r)  -  4/ (a;  1-  w)  +f  (x). 

But  the  given  function  is  of  the  third  degree.  Hence  for  it  we  have 
'^■'/(^')=0.     Noting  that  /'•  is  unity  we  obtain 

/(a:  +  4)  -4/(x  +  3)  +  6/(x  +  2)-4/(x+l)+/(a:)  =  0 
or,  denoting  /(a,'  +  n)  by  /„, 

/4-4/3  +  6/,-4/i+/o  =  0. 
Similarly  /^  _  4/, +  6/3-4/0  4-/1 -0. 
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Substituting  the  given  values  in  these  two  equations,  we  get 
/,  +  6/.  =  4-0809 
/,+   /,=  11819 
whence  /j  =  0-5798 ,  /,  =  0-6021, 

i.e.,         /(cc  +  2)  =  0-5798,      /(x-i- 4)  =  0-6021. 

Example  3. — If  when  a;  =  0,  1,  3,  4,  6,  the  corresponding  values  of 
/(x)  are  1,  1,  79,  253,  1291,  find  approximate  values  for  J'{x)  when  x  =  2 
and  when  x  =  b.     Under  what  conditions  will  the  result  be  accurate  ? 


3.    Expansion  of/(x  +  nw)  for  Integral  Values  of  n. 

We  shall  now  establish  certain  foiniulae  which  express 
f{x-\-mv)  (where  n  is  an  integer)  in  terms  of  differences  of 
various  orders.  These  formulae  will  be  shown  in  Chapter  II.  to 
be  valid  under  certain  conditions  even  when  >l  is  not  an  integer. 

It  will  be  found  convenient  to  introduce  another  operator  E 
defined  by  the  equation 

E=l+^. 

This  new  operator  will  obviously  obey  the  same  laws  as  A. 

Now     EJ  (x)  =  (1  +  ^)/ix)  =/{x)  +  A/(x)  =f(x  +  w). 

Thus  the  effect  of  the  operator  E  on  the  function  f{x)  is  to  increase 
its  argument  x  by  the  given  constant  quantity  w.  Repeating  the 
process  we  get 

E"J  (x)=f{x  +  n  w) 

and  since  the  operator  A  may  under  the  conditions  of  >:;  2  be  treated 
as  an  algebraical  quantity,  this  last  equation  may  be  written 

f{x  +  nw)^{\+^Yf{x) 

=/(x)  +  ,.c,  Ay  (x)  +  „c,  AV(x)  +   ...   +^'\f{x)  (1) 

where  „<\  denotes  the  coefficient  of  a''  in  the  binomial  expansion  of 
( 1  +  a)".  The  function  f{x  +  n  w)  is  here  expressed  in  terms  of 
/(•*-')  and  its  successive  differences. 

It  is  evident  that,  since  the  ef!ect  of  the  operator  E  on  f{x)  is 
to  increase  x  by  w,  the  result  of  operating  on  f{x)  with  E~^,  where 
E~^  is  defined  by  the  equation  EE~^  =  \,  will  be  to  subtract  w 
from  X, 

i.e.,  E~'/{x)=/{x-tv). 
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Hence  the  above  formula  (1)  may  be  transformed  as  follows  : 
y  (a;  +  nw)  =/{x)  +  „c,  A  fix)  +  (l  +  A)  E-'  {,.c„  A\f{x)  +  „c,  A"/(x) 

=/(.«)  +  „ci  A./  (x)  +  (1  +  A)  1,0,  Ar  f{x  -  w) 

+  „C3  AV(;^;  -  7/;)  +    ...    ^^'f{x-w)\ 

=J  ip)  +  nCi  Af{x)  +  „c.,  A-/(a:  -  w)  +  „+iC3  A"/ (■''•  -  w) 

+  (l+A)^-M,.+iC,Ay(x-M;)+   ...  } 

=/(^)  +  „Ci  A/ (a;)  +  „c.,  AVX*-'  -  i«)  +  -,+1^3  -^V(«  "  ^<^) 
+  (l+A){,.+iC,Ay(a;-2«;)+   ...  } 

=/(x-)  +  „Ci  A/(a;)  +  ,fi.  A"f{x  -  w)  +  „+iC,  Ay  (,x-  -  t«) 
+  „+,C4  A\/'(a3  -  2m;)  +  „+„c,  A'/(x  -  2m;)  +   ... 

Operating  next  with  (1  +  A)^-^  on  the  7th  and  subsequent 
terms,  then  on  those  beginning  with  the  9th  term,  and  so  on,  we 
obtain 

/{x  +  nw)  =f{x)  +  „Ci  Af{x)  +,,c,  A-f{x  -w)  +  „+,€.  Arf{x  -  tv) 

+  „^,c,  A'/{x  -  -Iw)  +  „^^c,  ^/{x  -2w)+  ...       (2) 

the  formula  ending  at  the  term  A'-"~^/(.x  -  n  -  I  iv). 

In  the  same  way  it  may  be  shown  that 

J\x  +  n  w)  =f{x)  +  „c,  Af{x  -w)  +  „+iCo  A-f{x  -  10) 

+  „^ic-,  Av(^  -  ^«^)  +  »+2C,  AV'{x  -  2*^-) 
+  „^,c-,,AV(x-- 3?,(;)+...  (3) 

the  formula  ending  at  the  term  A""_/'(;f  -  n  w). 

Changing  x  into  x  +  w  and  it  into  n  -  1  in  (3)  we  have 
f{x  +  n  w)  =/{x  +  w)  +  „_iCi  A/(a;)  +  ,.c.,  A-/(a;)  +  ,,03  AV"(x  -  w) 

+  n+ic,  AY(x  -w)  +  „+,c,  A'/{x  -2w)+...         (4) 
the  formula  ending  at  the  term  A-"J\x  -  11  -  2  w). 
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Adding  (2)  and  (3)  we  obtain 

A  f(x)  +  A  /"(.f  -  w)       n 
f{x  +  71  to)  =f{x)^  ,A  -^^^ :^ +  —  ,fix  A-/  {^  -  ^) 


A^/ix-w)  +  A'/(x-2w) 
■  2 

^^c,^\/\x-■2w)  +  ...  (5). 


Similarly,  addition  of  (2)  and  (4)  gives 

+  ..c,  ^'^''"^  +  -^V>  -  ">  ^  "^ .,,  ^VX.. ,  ^,  + . . .   (0). 

Formulae  (1),  (2),  (3),  (4),  (5),  (G)  may  be  compared  with  the  expansions 
introduced  in  Chapter  II.  under  the  names  Newton's,  Stirling's,  etc., 
Formulae  of  Interpolation. 

Again,  from  the  definition  of  the  operator  E,  we  see  that 

y'/{x)  =  (E-l)"/{x). 
Expanding  the  operator  in  this  last  relation  we  obtain 
AV(x)  =  Ey{x,  -  „c,  E-'Jl'-)  +  nC:  E'-"-J\x)  ...   +  (  -  1)V'{-*'') 
=J\x  +  n  w)  -  ,,cj\x  +  n-\w)  +  ,fi-if{^  +  r4  -  2  w) 
+    ...    +(-l/'/(x). 
In  this  result  we  have  the  ?tth  difference  of  the  function  f{x) 
expressed  explicitly  in  terms  of  /(x)  and  its  successive  values. 

4.    Relations   between   the    Dififerences    and    the    Dif- 
ferential Coefficients  of  a  Function. 

The  second  example  of  ^5  1   may  ha\e  reminded  the  reader  of  a 
somewhat  similar  result  in  the  differential  calculus,  viz., 

— -  sin  {ax->rb)  =  a"  sin  (ax  +  b  +  .1  /t tt). 
ax 

This  result  may  be  derived  from  the  example  in  question.      For  if 
we  put  e6'  =  A.i;,    where  Ax  denotes  a  small   increment    in  x,  and. 
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after  dividing  both  sides  of  the  equation  by  (A  x)",  make  A  x  tend 
to  zero,  we  get 

d"J  (x)  _   I      f2  sin  ^  a  .  A  X 


iX)         f      /-^  sin  *a  .  iAx    ■"   .      ,  ,       ,       ,        .  ^, 

x"  '—•I  Ax  ) 


=  a"  sin  (ax  +  b  +  ^nir). 

This  naturally  suggests  that  there  may  be  some  relations  connecting 
the  successive  differences  of  /(x)  and  the  successive  differential 
coefficients  of  the  same  function.  That  such  a  relation  does  exist 
we  now  proceed  to  show. 

Applying  Taylor's  Theorem*  to  the  function /'(«  +  ?t')  we  obtain 

fix  +  70)  =/(«)  +  tvf  (X)  +    '^/"  {X)+     ...     +    ^/*""  (^■)  +  •  ■  • 

or,  since  f{x  +  tv)  -f(x)  =  Af{x), 

Af{x)  =  tvf{x)+  |^,/"(«^)+   ■••   +  ^/'""  (*•)+••• 

which  is  a  relation  of  the  kind  sought  for. 
Again,  by  Taylor's  Tlieorem, 

j\x  +  tv)=f{x)+    tvf{x)+     ^/"(^■)+    ...    +     ^/<™' (,«)+... 

f{cc  +  -2iv)=f{x)  +  2wf{x)+''^^f"{x)+..    +li^'/<"')  («;)  +  ... 
and  generally 
fix  +  n  IV)  =f  {,■)  +  n  wf  (x)  +  ^^"  /"  {x)  +  ...+  ^^"  y  <'«'  (x)  +  ... 


Hence 


A"/{x)  =j\x  +  n  IV)  -  „Cj  f{x  +  n  -  1  iv)  +  „Cn/{x  +  n-2  tv) 

-    ...   +i-\)V(x) 

tv'"' 
=  2  — -  {n'"  -  ,.Ci  (71  -  1)'"  +  „Co  (n  -  2)'"  -   . . .  }./"<""  (x). 


*  It  would  be  inopportune  at  this  stage  to  introduce  a  discussion  of  the 
questions  of  convergence  which  arise  in  connection  with  the  application  of 
Taylor's  Theorem.  We  shall  assume  a  convergence  sufficient  for  our 
purpose. 
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But 

?i"'  -  „'•,  {n  -  1 )"'  +  „c,  (?i  -  2)'"  +   . . . 

=  0       if  ?u</t* 

=  nl      if  ni  =  n 

=  i  ?t  (?i+  1 )  !        if  m  =  ?t  +  1 

=  — •  ?t  (3  ?i  +  1 )  (?i  +  - )  !       if  m  =  71  +  2 
etc. 
Making-  use  of  these  results  in  the  expression  for  ^"./'{''-')  '^^e  find  that 

15  n-{n  +  1) 


6! 


^/;"+y "+"'  (x)  ■ 


Corollary  1.— Since  the  nth  differential  coefficient  of  a  function 
/(x)  of  the  7ith  degree  is  a  constant,  and  the  derivatives  of  higher 
order  are  zero,  we  see  again  from  this  expi-ession  for  A"/(->')  that 
the  «th  difference  of  a  polynomial  of  degree  n  is  a  constant,  and  the 
differences  of  higher  ordei-  zero. 

Corollary  2. — By  reversing  the  series  for  A"/(x-),  we  obtain  an 
expansion  for  the  «th  differential  coefficient  of  /(x)  in  terms  of  its 
differences,  viz. 

Corollary  3. — If  the  /itli  dillerence.s  of  ;i  function  are  constant, 
then  the  function  is  a  polynomial  of  the  nth  degree. 

For  since  the  wth  differences  are  constant,  the  higher  ditlerences 
are  zero,  and  therefore  the  ?tth  flifferential  coefficient  is  constant, 


*  Chrystal'B  Alfjebra,  Fart  ii.,  Chap.  xx\  u.,  g  •). 
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and   the  differential   coefficients  of  higher  order  zero.     Thus  the 
function,  when  expanded  by  Taylor's  Theorem,  takes  the  form 

/u:)  =f{X)  +  (X-  -  X)f  {X)  +  '"  "  f^'f"  (X)  +  ... 

n  ! 
and  this  is  of  the  n\h.  degree  in  x. 

Corollarii  4. — Since  e''=l+r—  +  -:r-:4-    ...    4 ;  +  ... 

•^  1  !      2  !  n ! 

we  see  that  the  series  which  repi'esents  f{x  +  w),  viz. 

f{x)  +  wf{x)+  —f"{x)  +  ... 
r ,  d      w-  /  d\-  iif  /  dY  1  ^  ,   , 

may  be  symbolically  repi'esented  as 

e"^/(x). 

Hence  Ef{x)  =  e'''^/{x). 

MISCELLANEOUS  EXAMPLES. 

(Assume  in  the  following  exercises  that  w  is  unity,  unless  otherwise  indicated). 

1.  Find  the  first  differences  of  the  functions 

tan  ax;    log  a: ;     (  a'-^  -  a*  )  /  (a  -  1 ). 

2.  Prove  that 

.  .  „  sin  ^ 


cosx  6  cos(a:+  1)^  ' 

3.  Prove  that 

1         _     I  d 

2- tan— 

4.  Find  the  ?*"'  differences  of  the  functions 

a;+l 
— — — — - .     a;m+n  ;     sinaa;coste. 

5.  Find  the  nth  difference  of  each  of  the  factorials  x(x  +  a)(x  +  2a)...  and 
there  being  m  factors  in  each  factorial  and  the  increment 


x{x  +  a){x  +  2a). 
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of  X  being  a:    and  show  how  to  develop  any  polynomial  in  x  in  a  seiies  f)f 
factorials. 

In  particular  show  that 

x'  =x  +  'lx(x-l)  +  6x(X'  l){x-2)  +  x{x  -l)(x-2)(a;-3). 

6.  If  /{x)  =  ay  +  i2'  prove  that 

A-'f(x)-3S/(x)  +  2f(x)  =  0. 

7.  Prove  that 

Ax"  ~  2\-x"  +  3  A^'x"  -  ={x-\)"  -   {x-2)". 

8.  Establish  the  formula  for  tl  e  9ith  ditierence  of  the  product  of  two 
functions 

A"  {/uf) .  4>  {.o  I  =  A"/(x)  .  0(.<)  +  n  d"-i/(.r  +  1) .  A  0  ( .,) 

?l(7l-l) 

+  \7-^"-V(^-  +  2) .  ^'<P{x)  +... 

9.  Show  that 

and  find  the  value  of 

A-/(.r)  +  —  (a  -  cos  b)  A  fix)  +  (l-  —  cosh+  -\)/(x) 
n  \         a  a-  / 

cos  b  X 
where  /{.r)  =—^^ — . 

10.  The  record  of  exports  for  the  year  1813  was  destroyed  by  fire.  Make 
an  estimate  of  their  value  given  that  the  values  of  ihe  exports  for  the  years 
1810,  1811,  1812,  1814,  1815,  1816  were  48,  33,  42,  45,  52,  and  42  million 
pounds  respectively. 

11.  It  is  asserted  that  a  quantity,  which  varies  from  day  to  day,  is  a 
rational  and  integral  function  of  the  day  of  the  month,  of  less  than  the  fifth 
degree,  and  that  its  values  on  the  first  seven  days  of  the  month  are  30,  30, 
28,  25,  22,  20,  20.  P]xamine  whether  these  assertions  are  consistent.  If  so, 
assume  them  to  be  true,  and  find  (1)  the  degree  of  the  function,  (2)  its  value 
on  the  sixteenth  day  of  the  month. 


[CH.  II 


CHAPTER     II. 

FORMULAE    OF    INTERPOLATION. 

5.    Introduction. 

The  problem  of  determining  the  value  of  a  function  for 
some  intermediate  value  of  the  independent  variable  or  argu- 
ment, when  the  value  of  the  function  is  given  for  a  set  of  discrete 
values  of  the  argument,  is  called  Interjyolation.  When  the  function 
obeys  some  definite  law  the  approximation  to  the  true  value  may  be 
made  as  accurate  as  may  be  desired  :  but  if,  as  often  happens, 
the  rigorous  analytical  formulae  which  occur  are  very  com- 
plicated, they  will  be  of  little  value  from  the  point  of  view  of 
computation.  It  will  then  be  necessary  to  use  a  method  analogous 
to  that  which  has  to  be  employed  when  the  properties  of  the 
function  are  entirely  unknown,  save  for  the  values  (usually  obtained 
as  a  result  of  experiment  or  observation)  which  it  has  for  a  certain 
limited  number  of  values  of  the  argument.  It  is  to  this  case  that 
our  attention  will  be  directed. 

When  the  observations  are  given  for  values  of  the  argument 
equally  distant  from  each  other,  and  are  in  arithmetical  progression 
so  that  they  can  be  plotted  on  a  straight  line,  the  possibility  of 
interpolating  is  obvious.  But,  in  general,  the  observed  values 
do  not  follow  any  apparent  law,  and  unless  some  assumptions 
are  made  regarding  the  form  of  the  function,  the  problem  is 
indeterminate.  It  is,  however,  customary  to  obtain  the  entry  for 
values  of  the  argument  not  far  distant  from  each  other  ;  and  when 
this  is  done,  in  most  practical  problems  it  is  found  that  all  the 
differences  at  some  particular  order  become  sensibly  zero.  In  such 
cases  the  observations  may  be  represented  with  considerable 
accuracy  by  a  polynomial  in  x  of  the  form 

y  =  flo  +  «!  a^  +  «-  *"  +    •  •  •    +  «n  *"• 
This  is  the  equation  of  a  parabola  of   the   nth.    order,    on    which 
account   the   method   now   to   be  described    is    sometimes    called 
Parabolic  Interpolation. 
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In   some  cases   it   is  advisable  to   calculate    the    function    not 
directly,    but    in    combination    with    some    other    function.       The 

calculation    of    the    values   of    the  probability   integral  e~^*-'' dt 

furnishes  an  example  of  this  kind,  since,  when  the  argument  becomes 
large,  the  integral  and  its  differences  become  inconveniently  small 
and  irregular.     In  this  case  it  will  be  found  convenient  to  compute 

the  function  e^-*-'  I    e~-^''dt.     At  other  times  some  function  of  the 


^-•^  f   e-i^''dt 


given  function  is  calculated ;  for  example,  if  the  observations 
formed  a  geometrical  progression,  then  we  should  prefer  to  make 
the  interpolation  on  their  logarithms  ratlier  than  on  the  quantities 
themselves. 

Wallis  (1616-1703)  may  be  said  to  have  laid  down  the  principle  of  imer- 
polation  in  his  approximation  to  the  area  of  the  circle  y  =  (x-  .c-)'  '.     In  this 

work  he  first  assumed  that  the  value  of  the  integral    I     (x-x'^f  '  dx  might  be 

Jo 
taken  as  the  geometric  mean  of  the  values  of 

I    (.r  -  .t-)"  cZ.r  and   I     {x-x-)dx,  i.e.  of  1  and  J. 

Subsequently  he  saw  that  this  was  not  exactly  true,  and  that  the  value  of 

I    (x  -  .v^y-  dx  must  obey  the  law  expressed  by  the  series 


i''--"'""-  l/^-"^'"'- [,' 


{x-x^dx; 


this  was  equivalent  to  interpolating  the  value  of   the  integral   under  con- 
sideration. 

6.    Lagrange's  Formula  of  Interpolation. 

We  shall  first  show  how  to  find  a  polynomial  of  degree  {it  -  \) 
at  most,  which  has  the  same  values  as  some  arbitrary  function  ./"(.r) 
foi'  n  distinct  values  of  x,  .say  x  =  a^,  a.,,  ...  a„.  This  polynomial 
will  be  represented  by  the  parabola  of  order  (?i  -  1)  already  men- 
tioned, and  may  be  expressed  in  the  form 

/(x)  =  Cq  +  c^x  +  c,  X-+  ...  +  c„_i «"-'. 
But  it   will  be  found  more  convenient   to   write   it,  as  is  always 
possible,  in  the  form 

fix)  =     i     ^r  («  -  «.)    •••    (•«  -  «>-l)  (•«  -  «.+l)    •••   (■'-  -  «n)  (1) 
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When  .'■  is  put  equal  to  a,,  all  the  terms  on  the  right-hand  side  of 
(1)  vanish  except  the  term  containing  J,.:  at  the  same  time /(a?) 
becomes  /'(«,)•     Hence 

A,  =/  (a^)  /  (a,.  -  ffli)  . . .  (a^  -  «,_,)  {a,  -  a^+0  . . .  (a^  -  a„) 
and  equation  (1)  becomes 

'^^■^^-,fl(a,-aj...(a,-a,_,)(«,.-«,.4i). ••(«>-«„)•       "^  ^ 

This  formula  expresses  f{x)  for  all  values  of  x  in  terms  of  its 
values  at  the  n  points  a^,  a.,,  . . .  «„.  If  /(x)  is  a  polynomial  of 
degree  less  than  n,  the  formula  gives  its  accurate  expression.  But 
if /(ic)  is  not  a  polynomial  of  degree  less  than  n,  then  it  is  no  longer 
possible  (in  default  of  other  information)  to  determine  f(x)  from  a 
knowledge  of  its  values  at  the  points  a^,  a„,  ...  «„ ;  for  if  any  one 
such  function  were  known  we  could  obtain  another  by  adding  to  it 
any  function  which  is  zero  at  the  points  a^,  a„,  ...  a,„  e.g.,  we  could 
add  the  function  {x  -  a^)  (,t  -  a.,)  ...  (x  -  a„)  ({>  (j-)  where  <^  [x)  does 
not  become  infinite  when  x  =  a-^,  a^,  ...  «„.  If,  however, /"(aj)  is  a 
function  of  which  nothing  is  known  except  that  it  takes  certain 
definite  values  corresponding  to  certain  definite  values  «],  a.,,  ...  «„ 
of  the  variable,  then  the  simplest  assumption  which  we  can  make 
regarding  its  form  is,  that  it  is  the  function  of  degree  (n  -  I)  which 
satisfies  the  given  conditions ;  and  this  function  is  given  by 
equation  (2).  The  expression  on  the  right-hand  side  of  (2)  may 
therefore  be  taken  as  the  expression  for  /(x)  for  all  values  of  x  in 
the  interval  considered.      It  may  be  written  in  the  form 

.,    s        "  ±(oc) ,.,     - 

r=i{x-a,.)(fi  (a,) 

where  <^  (x)  =  {x  -  a^)  {x  ~  a.2)  ...  (x  -  a,,). 

This  is  known  as  Lagrange's  Formula  of  Interpolation.  The 
advantage  which  it  has  over  other  formulae  for  the  same  purpose 
is  that  it  is  in  a  form  suitable  for  logarithmic  computation. 

This  latter  form  of  Lagrange's  formula  may  be  obtained  directly  from  the 
formula  previously  given.  It  can  also  be  obtained  very  readily  by  the 
Method  of  Partial  Fractions. 

For  consider  the  proper  fractional  function   ^  .   ,    where 

0  (X) 

<t>  [x)  =  (x-a{){x- ao)  ...{x-a,, ). 
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Decomposing  it  into  partial  fractions  we  have 

fix)  _    A,     ^     A^     ^  ^_  _^    An    ^ 
<p{.r)      .r  -  ci,     X  -  Uo      "     x-a,,' 

But,  by  the  general  theory  of  partial  fractions,  the  coefficient  A,-  is  given  by 

/{a,.) 
^'-  =  .'(«.)• 
Hence  .  '  \ 

/(.r)_         /{a,)  ^         f{a.^  ^        ,  /( «»  )  ^        /^ 

<p(x)     (x-ai)0'(ai)       (X  -  ao)  ((>' (ao)     '"     {x- an)<p'(an)  ^         V 

which  gives  r^, ,      ^ 

Example    /.—Construct   a  polynomial   of   the   third    degree  in    .r   which   ^/ 
shall  have  the  values   -5,  1,  1,  7  when  the  values  of  x  are   -2,    -1,   1,  2 
respectively,  and  determine  its  value  when  x  is  f. 

Let  f{x)    be  the    required    function.     Substituting    the    given  values    in 
Lagrange's  Interpolation  Formula  we  have 

(x+l)(x-\)(x-2)  (x  +  2){x-\){x-2) 

Ji^'l  -(_2+l)(-2-I)(-2-2)^     ^'^  (-1+2)  (-1 -!)(-!  -2)^  ''^' 

(X  +  2)  {X  +  I)  {x  -  2)  {x  +  2)  {X  +  I)  {X  -   1) 

(1  +  2)  (1  -(-  1)  (1  -  2)  ^  (2  +  2)  {2  +  1)  (2  -  1) 

=  1%  (.r  +  1)  {x  -  1)  (x  -  2)  +  ^  (x +  2)  (X  -  1)  (X  -  2) 

-  h  (X  +  2)  (X  +  1 )  (X  -  2)  + 1.  (X  +  2)  (X  +  1 )  (X  -  1) 

•     =X3-X+I. 

From  this  it  follows  that/(§)  =  2|.  But  in  practice  it  is  not  necessary  to 
obtain  the  algebraic  form  of  f(x)  if  we  merely  require  its  value  for  some 
value  of  X.  Direct  substitution  in  the  formula  will  give  the  result  much 
more  quickly.     In  the  question  under  consideration 

1|+1)(-|-1)  (i-2)  (1  +  2)  (S-1)  (f-2) 

./(:j)-(_2+i)(_2-l)(-2-2)<     5)+    ^  _i  ^g)  (  -  1  -  I)  (  -  1  -  2) 

(5  +  2)  (^+1)  (g-2)         (§  +  2)  (g  +  1)  (§-1) 
+   (1+2)  (l  +  I)  (1-2)    ■*"  (2  +  2)  (2  +  1)  (2+1)    ^'^'' 

—  01 

Example  5. —Find  the  minimum   value  of  the  polynomial  of  the  second 
degree  which  has  the  following  values 

0  ,         1  ,         5 
/(x)       11,         5,       21. 

G.  2 
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Substituting  these  values  in  Lagrange's  Interpolation  Formula  we  have 
(x-l){x-5)  X  (x  -  5)  x{x-\) 

f  ir)  =  -TV.-^rr  -11  +  T— :T  •  ^  +     5T4"-  •  ^i 

.-.        f(x)  =V(2:f-6)  -  i(2.T-5)  +  f^(2.T-l). 

For  a  minimum  turning  %-alue  /'  (a;)  =  0 

i.e.,         44  (2  .r  -  6)-  25  (2  .r-5)  +  21  (2  .t-  1)  =  0 
80.T- 160  =  0 


Hence  the  minimum  value  of /(.r)  corresponds  to  the  value  2  of  the  argument. 
The  corresponding  value  of  f{x)  is  ,'?. 

Example  3.—lif{x)  is  a  function  of  .i*  which  assumes  the  values  20,  .35,  56,  84, 
when  X  has  the  values  0,  I,  2,  3  rospectivel}',  express  /(:r)  in  terms  of  x,  on 
the  supposition  that  M^f(x)-0.  Also  find  an  approximate  value  of  x  if 
/(.r)  =  25,  third  differences  being  neglected. 

Example  4. — If  four  consecutive  values  of  a  function  corresponding  to  the 
values  1,  2,  3,  4  of  the  variable  be  601,  777,  902,  999  respectively,  find  the 
value  of  the  fi;nction  when  the  argument  has  the  value  1'5. 

7.  Periodic  Interpolation. 

When  there  is  reason  to  suppose  that  the  function  is  periodic,  it  is  better 
to  assume  lor  f{x)  an  expansion  in  terms  of  periodic  functions,  say, 

f(x)  —  a(j  +  a■^  cos  a:  +  ag  cos  2  ic  +  ...  +  a„  cos  n  x 
+  61  sin  X  +  63  sin  2  .r  +  ...  +  6„  sin  n  x, 

the  number  of  unknown  coefficients  corresponding  to  the  number  of  observa- 
tions given.  When  the  observations  are  given  for  values  of  the  argument 
spaced  at  equal  intervals,  the  problem  becemes  one  in  Fourier's  analysis  and 
will  be  found  treated  in  the  tract  on  that  subject.  But  when  the  values  of  x 
are  not  in  arithmetical  progression  we  may  use  a  formula  which  is  analogous 
to  Lagrange's  formula  of  interpolation,  and  which  Gauss  has  shown  to  be 
equivalent  to  the  above  expansion  for/(.T). 

If/(a;)  is  given  for  x—Xq,  x-^,  Xo,  ...  a,'2„  the  corresponding  formula  is 

ff  \      ^    sm^(x-XQ)  ...sin^{x-Xr-^)s\n^{x-Xr+i)...am^{x-X2n)     ^,     , 

/W  =  2j f{Xr). 

r-O  Sin  J  {Xr  -  .Tj)  . . .  sin  ^  (.T,.  -  .Tr-l)  siu  ^  {Xr  -  Xr+l)  . . .  sin  ^  (Xr  -  X^n) 

8.  Notation. 

Hitherto  the  points  a^,  a^,  ...  a„  have  been  taken  to  be  any 
points  whatever.  But  in  the  most  important  class  of  cases  it  is 
possible  to  obtain  the  entry  for  values  of  the  argument  equally 
distant  from  each  other.     In  what  follows  we  shall  assume  that  the 
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ur^'uments  are  spaced  at  equal  intervals,  and  shall  use  the  notation 

of  tlie  foUowini;  scheme  which    for  many  purposes  *  is  more  con- 
venient than  that   used  in  Chapter  I. 


Argument. 

Entry. 

A 

A- 

A^                 A* 

A= 

rt  -  3  )(• 

^A-i 

a-2w 

y'-2 

i^^A--2 

SV- 

2 

/^./■-= 

S/-5 

H-^V- 

^           3'/-^ 

a  -  ir 

./■-I 

h¥-i 

«V- 

1     /x^V-i       3V"-i 

l\f-l 

^A     i-^Vi 

sy_ 

a 

A 

/-sv; 

Rt\ 

sy^ 

a  +  ir 

/. 

5/J 

\>-h 

a-\-'2  IV 

A 

f^^A 

8v; 

H-A 

a  +  Zw 

A 

when* 

¥-h 

=./;-/- 

1 

f^Ai  =  h  (Ao  +  /i) 

¥h 

=./;-/n 

etc. 

etc. 

^) 

sVo 


9.    Newton's  Formula  of  Interpolation. 

If  we  assume  that  the  differences  a.,  -  a„  a.^-  a^,  ...  in  Lagrange's 
int(>rpolation  formula  are  all  equal  to  tr,  and  that  a,  in  the  formula 
corresponds  to  a  in  the  scheme,  then  the  quantities  /(a,.)  take  the 
following  forms  :  — 

/■(«3)  -A  =y; + 8/?  =  /, + 2  8/j  +  5v. 


/( 


a.o=./;-^/..^^^^v;.^^^^::^^^^ 


Especially  in  connexion  with  central-difference  formulae. 
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Lagrange's  Tntcrpolation  Formula  then  becomes,  if  the  number 
of  given  points  is  denoted  by  k, 

_     {x  -  «.)  (a?  -  as)  {x  -a,)  ...  (x-  %) 
•^  ^""^  "  {a,  -  a,)  (a,  -  a,)  {a,  -a,)...  (a,  -  a,)  ^' 

{x-a,){x-  g,)  (x-a,)  ...{x-  a,)  +8/\) 

{a„  -  «i)  {a.y  -  a;j)  (a^  -  ^j)  •  •  •  («2  -  «t) 
.      (.V  -  a,)  (.r  -  a,)  (x-  a,)...  {x  -  a,)    ^^^.  ^  g  g^-    ^  g.^^ 


(a,  -  a,)  {as  -  a..)  {a,  -a,)  ...  {gs  -  «,) 

(x  -  ttj)  . . .  (a;  -  a,)  (x  -  a^^.)  ...  («  -  a^,) 


Avhere  a^  =  a,  a.,  =  a  +  w,  . . .,  a,,  =  a  +  (^  -  1 )  ?^. 

This  is  an  expression  of  degree  (^  -  1 )  in  x.  In  it  the  coefficient 
of  y^  is  unity  when  x  =  a„  a^,  ...  a^,  and  hence  is  unity  always. 
Again,  the  coefficient  of  Sy^  is  0  at  a^,  1  at  a^,  2  at  ag,  ...  r  at 
a,._^i,  ...,  and  as  this  coefficient  must  be  a  polynomial  in  x  of  degree 

less  than  k,  it  must  be  no  other  than .      In  fact  the  terms  in 

w 

8/\    in    the    above    series    simply   constitute    the    expansion    of 

8/j   by  Lagrange's  formula.     We  shall  denote  by  n,  so 

w        ■^'  w 

that    n    is    a    real    quantity,    not   necessarily   an    integer.       The 

coefficient  of  8/^  is  thus  n.     Similarly  the  coefficient  of  S'-yj  is  0 

r(r-l) 
at  «!,  0  at  ao,  1  at  a^j,  . . .  — - — —  at  «,.+,  . . . ,  and  in  the  same  way 

we  see   that  this   coefficient  is      -^^^-j — -  for  any  value  of  n.     The 

other  coefficients  can  be   determined    in    the   same    way,    and    we 
thus  finally  have  Lagrange's  formula  transformed  into  the  equation 
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This  is  called  Newton^s  Formula  of  Interpolation*  In  it  the 
differences  are  those  which  occur  along  a  line  parallel  to  the  upper 
side  of  the  triangle  in  the  scheme.  It  may  on  this  account  be 
termed  a  Forward  Difference  formula  of  interpolation,  and  will  be 
found  very  useful  when  the  value  of  the  function  near  the  beginning 
of  a  limited  series  of  observations  is  required. 

It  must  be  understood  that  the  function  furnished  by  Newton's 
formula  is  the  polynomial  of  lowest  degree  which  has  the  given 
values  at  the  k  given  points  a,  a  +  ic,  a-{-^w,  ...  ,  a  +  {k-\)w: 
and  the  series  on  the  right  is  a  terminating  series.  If  the  given 
values  at  the  k  given  points  are  the  values  of  some  given  function 
4>  {x)  at  the  points,  we  cannot  conclude  that  <^  {x)  is  represented  by 
Newton's  formula  unless  4>  (*)  is  the  polynomial  of  lowest  degree 
which  passes  through  the  points. 

Suppose,  now,  that  the  values  of  a  function  are  given  not 
merely  for  a  finite  set  of  values  of  its  argument,  but  for  the  infinite 
set 

,>;  =  a-\-r  IV  (r  ^  U,  1,  2,  3,  . . .  ao?  inf.). 

Then  we  can,  by  Newton's  formula,  construct  a  polynomial  of 
degree  (k  -  \)  which  has  the  given  values  for  the  first  k  of  these 
values  of  the  argument,  where  k  is  any  positive  integer :  and  it 
may  happen  that,  as  k  increases  indefinitely,  the  series  on  the  right 
in  Newton's  formula  becomes  (for  some  range  of  values  of  x)  a  con- 
vergent series,  having  for  its  sum  a  function  f{x).  The  function 
f{x)  will  now  have  the  given  values  for  each  of  the  values  a  +  rw 
(r  =  0,  1,  2,  3,  ...  ad  inf.)  of  the  argument.  There  are,  however,  an 
infinite  number  of  functions  which  satisfy  this  condition  and  f{x) 
is  that  one  of  them  which  is  the  limit,  when  A;  — >  oo  ,  of  the  polynomial 
of  lowest  degree  which  has  k  of  the  given  values. 

An  expression  for  the  difference  between  the  value  of  the 
function  and  the  sum  of  the  first  (r+1)  term.s  of  the  limiting 
form  of  this  j)olvnonTial  will  now  be  obtained.      Let 

n(ii-  \)  ...  (n-  r  +1)        . 
/  (a  +  u  IV)  =yo  +  "  <^yA  +  •  •  •   +  —^^ (>Jj:   +  li- 

Then 

m{ni  -  \)  ...  {m  ~  r -\-\)  ^     . 
f{a  +  m  w)  -Jq  -  m  S /^  -   . . . j 6',/^    -  K 


*  Cf.  §  3,  Formula  (1). 


INTERPOLATION  AND  NUMERICAL  INTEGRATION      [cii. 


will  be  zero  for  m  =  0,  1,  2,  ...  r  and  n  :  its  derivate  of  order  (r+  1) 
will  therefore  be  sero  for  some  value  7i  of  m  lying  between  the 
least  and  the  greatest  of  the  quantities  0,  r,  ?i,  and  hence  the 
residue  is  given  by 

It  may  be  remarked  that  Newton  discovered  the  binomial  coethcients 
when  studying  the  expansion  in  series  of  a  function  by  interpolations.  He 
first  of  all  considered  the  expansions  of  the  expressions  (1 -:c-)°,  (l-.i'-)\ 
(I  -x^f,  ...  and  deduced  from  them  the  expansions  of  (l  —  x"^}^,  (1  -x^)^ ,  .... 
Then,  by  analogy,  he  obtained  the  expression  for  the  general  term  in  the 
expansion  of  a  binomial,  and  thus  established  the  binomial  theorem. 

Example  1. — From  the  following  table  of  values  of  the  function  /(.v) 
determine  the  value  of /(3-5). 


X                 fix)                      A                     A--' 

A3 

1             0-208460 

29242 

2             0 '237702                                  -213 

29029 

27 

3             0-266731                                  -240 

28789 

■26 

4             0-295520                                  -266 

28523     - 

5             0-324043 

In  this  example                    a-1,     n-^,     >r-l 

/„-/(!)  =  0  208460. 

Hence                        f(i-o)=Jl  f  ^5/^  -  15-/1  +  A  oVs. 

=  0-208460 

14621 

27-2 

=  0-2-23106. 

Example  2. — Corresponding  values  of  x  and/{,t)  are  given 

in 

table  :— 

fix) 

4  0                54-5982 

4-1                60-3403 

4-2                66-6863 

4-3                73-6998 

4-4                81-4509 

4-5                90-0171 

Obtain  the  value  of /(4-05). 
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A  certain  amount  of  discretion  is  necessary  in  the  application 
of  Newton's  formula.  Let  us  assume,  for  instance,  that  the 
given  observations  are  finite  —  say  seven  —  in  number,  and 
correspond  to  the  values  a-Sw,  a-2w,  ...  a  +  3iv  of  the 
argument  as  in  the  scheme,  §8.  Ji  a-\-nw  lies  between  a-3tv 
and  a-2w,  then  Newton's  formula  ior/{a  +  nw)  will  involve  the 
quantities. /I3,  8/_,,  S-/_„,  8y_ ,,  0^1],  ^/_y  and  S\f;„  all  of  which 
occur  in  tlie  scheme.  Thus  a  good  approximation  to  the  value  of 
/(a  +  mv)  may  be  obtained  by  means  of  the  formula.  But  if 
a  +  mo  lies  between  a  +  2w  and  a  +  3  w,  then  in  the  Newton's  formula 
iov  f{a  +  nw)  the  only  known  terms  are  the  first  two  depending 
on  J\  and  8/i;.  In  general,  these  are  quite  insufficient  to  determine 
the  value  of  /{a  +  n  w)  to  the  necessary  degree  of  accuracy.  It  is 
therefore  expedient  in  this  case  to  find  a  formula  involving  a 
larger  number  of  known  quantities.  A  suitable  one  may  be 
obtained  by  replacing  iv  in  the  original  formula  by  -  w.  Then, 
since 

we  must  replace  0 /•  by 

Similarly,  since 


we  must  replace  (>-J\  by  \ 

./o  -  2./L,  ^JL;  or  d-./L, . 

Tn  the  same  way  we  can  determine  the  functions  which  must 
take  the  places  of  8V|)  ^^.t\ Newton's  formula  then  becomes 

'/i  (n  -  1 )  ,.,   .        ?t  (n  -  1)  (m  -  2) -.  . 
J{a-nw)  =/,  -n8j_,+     \^  ^        crj_,  -  -^ ^^- S'V  -;;+••  ■ 

This  expansion  involves  differences  along  a  line  parallel  to  the 
lower  side  in  the  scheme,  and  it  may  therefore  be  regarded  as 
a  formula  for  Backward  Interpolation.  When  a  +  n  w  lies  between 
a  +  2w  and  a  +  3  w,  the  formula  will  depend  on  the  (juantities 
A,  Syj,  8-/;,  8V§,  SV;,  Syj,  and  SV„,  all  of  which  are  known  :  it 
is  therefore  very  convenient  when  the  value  of  /(x)  is  required  for 
values  of  .';  near  the  end  of  the  series  of  observations. 
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Example  5.— From  the  following  table  of   values  of /(x)  determine  the 
values  of /( 104 -25). 

A3 


.(■ 

/(•••) 

A 

A- 

101 

IU301'J8 

30906 

102 

1061104 

31518 

612 

103 

1092622 

32136 

618 

104 

1124758 

32760 

624 

105 


1157518 


In  this  example    a  -  105,     n  =  0-75,     *'•  =  1 

/„=/(105)  =  1157518 
/(104-25) 

0-75  (0-75-1)              0-75  (0-75-1)  (0-75 -2) 
=  /o-0-75§/_,+- oy_j- ~ 5y-5H 

=/o  -  0 -75  5  /•_ ,  -  0 -09375  5-/_j  -  0  -0390625  S^f.  +  . . . 
=  1157518-0-75  (32760)  -0-09375  (6-24)  -  0-0390625  (6) 
=  1157518-24570 
58-5 

0-234375 
=  1132889-265625. 
Example  4. — Using  the  data  of  Example  2,  find  the  value  of/(4-45). 


10.  Stirling's  Formula  of  Interpolation. 

As  Newton's  formula  involves  differences  which  lie  along  a 
slanting  line  it  will  frequently  not  be  suitable  for  calculating  the 
value  of  the  function  for  a  value  of  the  argument  near  the  middle 
of  a  series  of  observations.  By  a  suitable  transformation  of  this 
formula,  however,  we  may  derive  another  involving  differences 
which  lie  on  a  horizontal  line  ;  such  a  formula  is  termed  a  Central- 
Difference  Formula. 

Since 


we  have 
Hence 


8ji  =/xo/;  +  ^S7;. 
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In  the  same  wav  the  higher  diffeiences  which  occur  in  Newton's 
formula  may  be  determined  in  terms  of  the  diiferences  which  occur 
on  the  same  hoi'izontal  liueasy,",.  Substituting  these  in  Newton's 
formula  we  obtain 

J  {a  +  n  iv)  -^Jo  +  n[ix  d/^  +  l  ()-J„\  +  — { d-J^,  +  p.  6-J„  +  i  6Vo} 

71  (n  -  1 )  ill  -  2)  ,     , ^,   .         ,.   .      ,  OK  • . 

+  -^ — ^} ■  {/^  o'/o  +  ^  ^\r,  +  II  (V7„  +  ^-  ?>%]  +  ... 

«'  ,,   .       n(n'~\)      ,.,  .       n-(n--\)„,   . 
=Jo  +  n  1^  ojo  +  71  67o  +  3,        1^  o;/o  +  -^Jl 5Vo 

+  ""''-^>,<"^-^V8-/>  ... 

This  is  known  as  Stirling  s  Formula  of  Interpolation.  It  should 
not  be  used  for  determining  the  value  of  a  function  for  values  of  the 
argument  near  the  beginning  or  the  end  of  a  series  of  observations, 
as  it  would  not  then  be  possible  to  obtain  a  sufficiently  high  order 
of  differences. 

E.xample  1. — Using  the  data  of  §  9,  Example  1,  determine  the  value  of 
/(3-5). 


.<■ 

/U-) 

A 

A- 

A^ 

1 

0-208460 

29242 

2 

0-237702 

29029 

-213 

-27 

3 

0-266731 

(28909) 
28789 

-240 

(-27) 
-26 

4 

0-295520 

(28656) 
285-23 

-266 

(-26) 
-•26 

-) 

0-324043 

-  292 

•28-231 
6  0  •35-2274 

In  this  case     't  —  3,     "•  =  1,     n  —  h 
/„= /(3)  =  ©••266731 . 
/(3-5)  =/o  +  ^  M  5/0  H-  \  5V,  -  iV,  5-\/-o 
^  0  •266731 
14455 

2-30 
=  0-281158. 
Example  J.  — Using  the  data  of  §  9,  Example  2,  tind  the  value  of /(4-21). 
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If  n  is  greater  than  05,  it  is  better  to  use  backward  differences. 
The  corresponding  formula  may  be  obtained  by  merely  changing 
the  sign  of  w  in  Stirling's  Forward  Difference  Formula,  viz.  :— 

/{a  -  n  w)  = /o  -nil  6 /,  +  —  8-Jo ^-^ /'-  SVo 

Of  course,  either  formula  may  be  applied  whether  7i  is  greater  or 
less  than  0-5  :  thus,  when  a  check  is  necessary,  it  will  be  found 
better  to  compute/(a:)  by  the  formula  not  already  used. 

Example  5.— Use  this  Bickward  Difference  Formula  to  check  Example  L 
Here  a  =  4,     v  =1,     n  =  h 

/o  =/(*)  =  0-295520 
/(3-5)  =/u  -  4m5/o  +  i  S^-f,  +  iVSVo 
=  0-295520  -  14328 


=  0-281157. 

The  agreement  with   the   former   result   is   perfect  when   we  take   into 
consideration  the  fact  that  the  last  figure  is  forced  in  both  cases. 

Example  .^.— Using  the  data  of  §  9,  Example  2,  find  the  value  of/(4-29). 

11.    Gauss'  and  Bessel's  Formulae  of  Interpolation. 

Another  central  difference  formula  may  be  obtained  by  trans- 
forming Stirling's  formula. 

Since  fi  8jl  =  8/^  -  -|  S'-Jl 

/x6v;  =  8v;-i5v;       etc., 

we  see  that  Stirling's  formula  may  be  written  in  the  form 

J{a  +  n  w)  =J,  +  n  8j^  +      ^^,     '  SVj,  +  ^ '-jf '-  8\/l 

(n+l)7i{n-l){n-2) 

^  (n+2)(n+\)7i{7i-l){n-2)  g^^.  _^  ^ 

In  this  result,  which    is    frequently    known    as   Gauss'    Formula, 
the  /A-terms  do  not   occur.      The  differences  which  occur  in  this 
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formula  are  those  which  lie  on  two  adjacent  horizontal  lines  of 
the  scheme  in  §  8. 

Again,  since  t\  +/» =  -  IKt'^ 

we  have 

5Vo=/^5^./;,-isv;     etc., 

and  hence 
j\a  +  nw) 

_^  {n  +  2){n+l)n{u-l){7i-2)  ^^  ^.  _^ 

5  !  •  i      ■■■ 

,^-  .       n(n-l)      -,          n  (71  -  I)  (n  -  h)  r,.,  ,. 
=  fij,+{7i-h)8j^  +     ^^,     V^71+  -^ 3^r ^SV^ 

(-M+l)7^(n-l)(n-2) 

4  :  '     '  i 

(«+l)  u(m-  1)  (/^  -•-')  (h  -  .V)  ,.  , 
+ — —  6'/,  +  ... 

D  ! 

This  is  frequently  known  as  BesseUs  Formula.  In  it  the  differ- 
ences are  those  which  occur  on  a  horizontal  line  mid-way  between 
the  lines  on  which /„  and/i  stand. 

The  formula  most  suitable  for  the  construction  of  mathematical 
tables,  when  differences  of  higher  order  than  the  fourth  arc 
neglected,  is  obtained  from  the  first  of  these  on  putting 

When  this  substitution  is  made  we  get 

n{\  -n)  ^.   .       ?t  (I  -  n-)  ...  . 
J  (a  +  n  iv)  =yo  +  «  6 /^  +      \,  ,       o-y, ^yj —  6;/^ 

n(\  -  ri-)  (2  -  n)      ^.   . 
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Bessel's  formula  will  be  found  particularly  suitable  when  n  is 
equal  to  -^,  for  in  this  case  the  formula  becomes 

in  which  only  even  differences  occur. 

As  with  Stirling's  formula,  none  of  these  results  should  be 
employed  when  the  value  of  the  function  is  required  for  values 
of  the  argument  near  the  beginning  or  the  end  of  a  series  of 
observations.  Stirling's  formula  will  be  found  more  convenient 
when  71  is  very  small,  or  when  the  number  of  observations  given 
is  odd  :  Bessel's  when  n  is  nearly  equal  to  a  half,  or  the  number 
of  observations  even. 

The  backward  difference  formula  corresponding  to  Bessel's 
forward  difference  formula  may  be  obtained  in  the  same  way  as  in 
Newton's  and  Stirling's  cases.  But,  except  as  regards  sign,  the 
coefficients  in  the  two  formulae  will  be  the  same,  and  therefore  the 
one  cannot  be  used  as  a  check  on  the  other.  When  a  check  is 
required  to  a  result  computed  from  Bessel's  formula,  it  is  better 
to  compute  it  again  using  Stirling's  formula. 

Example  1. — The  values  of  a  function  for  the  values  -2,  —1,  0,  1,  li,  of 
the  argument  are  respectively  0-12569,  0-17882,  0-23004,  0  27974,  0-3-28-23. 
Using  Gauss'  formula,  show  that  when  the  argument  has  the  value  0-4,  the 
value  of  the  function  is  0-2500S. 


X 

f(x) 

A 

A- 

_.2 

0-12569 

5313 

-1 

0-17882 

5122 

-191 

0 

0-23004 

4970 

-152 

1 

0-27974 

4849 

-121 

2 

0-32823 

In 

this  example 

a  =  0,     /( 

=  0-4,     w  =  1 

31 


X,  =/(0)  =  0-23004. 
Hence  /(0-4)  =/«  +  0  4  S/.  -  0-12  S'-f,      0-056  5Vi 

=  0-23004 
1988 
18-2 
=  0-25008. 
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Example  2. — Using  the  data  of  §9,  Example  2,  find  the  values  of/(4-21) 
and /{4 -29). 

12.    Evaluation  of  the  Derivatives  of  the  Function. 

When,  as  often  happens,  the  values  of  the  differential  co- 
efficients of  a  function  which  is  given  by  a  series  of  observations 
are  required,  use  may  be  made  of  the  symbolic  equality  cstabHshed 
in  §  4,  Cor.  4. 

For,  since 

e''^\f{x)  =  Ef{x) 
we  have 

w^^.f{x)  =  \o2,E.f{x). 

The  method  will  be  evident  from  the  following  example  : — 

Example  1. — The  values  of  a  function  corresponding  to  the  values  TOl, 
1-02,  1-03,  104,  1  05  of  the  argument  are  1  030301,  1  061208,  1  092727, 
1 -124864,  1  •157625  respectively.  Find  the  value  of  the  first  diff'erential 
coefficient  of  the  function  when  the  argument  has  the  value  1  -03. 

Here  -'•  =  0-01.         Let /o  =/(  1-01). 

d 
Then  w  j^J{  1  -03)  =  log  £ .  /( 1  03) 

=  \ogE./o 
^  log  E.E'f,, 
=  log(l  +  A).(l  +  A)'-^../;, 
/  A2         A3  X 

Neglecting  difi'erences  of  higher  order  than  the  4th,  this  gives 
0-01-^/(1  -03)  =  { A  +  S  A2  +  i  A^'  -  ,V  A^)/„ 

=  5(/3-/i)-iV(/4-/o) 
=  0-042437 -0-010610 
=  0-03 18-27 
/'(r03)  =  3-1827. 

Example  .?.— The  values  of  a  function  corresponding  to  the  values  15,  2, 
2-5,  3,  3-5,  4,  of  the  argument  are  4-625,  2-000,  -0125,  -1,  0125,  4  000, 
respectively.  Find  the  value  of  the  first  differential  coefficient  of  the  function 
at  the  points  where  the  argument  has  the  values  2-5  and  3-0. 
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The  method  is  perfectly  general.  The  same  result  may  be 
obtained  by  differentiating  any  of  the  formulae  of  interpolation 
already  established.  For  example,  if  we  differentiate  with  respect 
to  n  the  formula  of  Stirling, 

A-j-  ?i  iTi" 1 ) 

/(a  +  n  w)  =./'o  +  n  1x8/1+  .77  ^'/o  +   3-j —  H-  ^"'/o 

,    n-(^r-l)  , 

o  Jqt   ... 


4  ! 
we  have 

w/'  (a  +  71  w)  =  /x  8/1  +  n  8- /l  + ^  {3n-  -  I)  fx  8% 
+  ^\{2'n;'-n)8%+  ... 
w-f"  {a  -t- n  iv)  =  S-'./o  +  nix 8\fo  +  i\  (6  rtr  -  1 )  8*/l  +  ... 

Putting    M  =  0  these  become 

fvf'  (a)  =  [JL  8/1  -  i /u, S\tl  +  ;/o  fi  8%  -  ... 
,v\r{a)  =  8-%-^\8%  +  ,}^8%-... 

These  formulae  are  of  use  in  determining  the  turning  values  of 
a  function,  the  values  of  the  argument  corresponding  to  given 
values  of  the  differential  coefficients,  the  points  of  inflexion  on  a 
curve  of  which  isolated  points  are  given,  and  also  in  determining 
whether  a  series  of  observations  is  periodic. 

Example  3. — The  consecutive  daily  observations  of  a  function  being 
0-099833,  0-208460,  0-314566,  0-416871,  0-514136,  0-605186,  0-688921, 
0*764329,  show  that  the  funclion  is  periodic  and  determine  its  period. 

From  the  given  observations  we  obtain  the  table  :  — 


X 

1 

0  099833 

A 
108627 

A- 

A-' 

A* 

2 

0-208460 

106106 

-2521 

-1280 

3 

0-314566 

102305 

-3801 

-1239 

41 

4 

0-416871 

97265 

-5040 

-1175 

64 

5 

0-514136 

91050 

-  6215 

-1100 

75 

6 

0-605186 

83735 

-7315 

-1012 

88 

7 

0-688921 

75408 

-8327 

8 

0-764329 
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Taking  a  =  3,    /o  =/(«)  =./"{3),     we  have,  since  w  =  \, 

=  -3801-ii5-41 
=  -3804 
1  3804 


•   •       /{3)  --^    ^^^  -       314566 

=  -  0  0121. 

In  the  same  way  we  find  that 

(i)     When           «  =  4,    /o=/(4) 

/"(4)  =   -5045,     -^f 

(4)  -  -0-0121 

(ii)     When         a  =  5,    /o  =/(5) 

/"(5)=-6-221,    ^^-:^/ 

'(5)=   -0  0121 

(iii)     When        a  =  6,    ,t\  =  /(6) 

/"(6)=  -7322,    ^-/ 

'(6)  =   -0-0121 

Hence  in  all  cases  j^^ /"  (.r)  =  -0-0121. 

Putting  y  —  J 

differential  equation 

1     d-y 

-0-0121 

f(.c),  this  gives  the 


whence  ^^,  +  0-0121  y  =  Q 

y  -  A  cos  0-11  .T +  5  sin  0-11. T. 
This  result  shows  that  y  or ./'(.')  is  a  periodic  function  of  .»',  and  that  its 
period  is  27r/011,  or  5712  days. 

Example    ^.— The    following    table    represents   one   of   the   Bessel   func- 
tions, /   (.r)  : — 


1-0 

V 
0-765198 

ri 

0-719622 

1-2 

0-671133 

1-3 

0-620086 

1-4 

0-566855 

1-5 

0-5118-2S 

1-6 

0-455402 

dy         dhf 
Find  the  values  of  -j^  and  ^„ 

for  .1 

=  1-3,  and 

four-place  accuracy)  in  the  equation 

dx-  ^    X 

dy 

(.-$) 

1-3,  and  hence,  by  substituting  (with 


determine  the  value  of  7i. 
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Forming  the  table  of  diifereuces  we  have 


X 

1-0 

2/ 
0-765198 

A 
-  45576 

A- 

A3 

11 

0-719622 

-  48489 

-2913 

355 

1-2 

0-671133 

-51047 

-2558 

374 

1-3 

0-6-20086    ( 

-52139) 
-53231 

-2184 

(381) 
388 

•1-4 

0-566855 

-  55027 

-1790 

397 

1-5 

0-511828 

-  564-26 

-1399 

1-6 

0-455402 

Here 

w  =  0-1, 

,  a  =  1-3, 

A=f(a) 

=  /(l-3). 

Hence         0-1/' (I -3)  =  -0-052139  -^0-000381) 
/'(1-3)  =  -0-52203 
•01/"(l-3)  =   -0-002184  -1^(0-000014) 
/"(1-3)  =  -0-2185. 
Substituting  these  in 

d-y        I    dy       /  n-  \ 

we  have 

-0-2185  +  j^( -0-5220)  +  (l  -  j^^  (0-6201)  =  0 

-0-2185  -  0-4016  +  (l  -  Jf^i)  (0-6-201)  =  0 

-0-6201  +  (l  -  lY^i)  (0-6201)  =  0 

11^  =  0 


1  +  1 


Hence  the  Bessel  function  in  question  is  J^ (a). 

Example  5. — Show  that  the  following  observations  are  these  of  a  periodic 
function,  and  find  its  period  : — 


X  (years) 

./•(a-) 

1 

0-198669 

2 

0-295520 

3 

0-389418 

4 

0-479425 

5 

0-564642 

6 

0-644217 

7 

0-717356 

8 

0-7833-27 
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Example  6.  -Observations  were  taken  of  tlie  cooling  of  a  boiler,  and  the 
results  are  given  in  the  following  table  :  - 

t     (time  in  hours)  1-2  19  •2-6  :V3  40 

e    (temperature   F)  182-3         175-2         ITl'l         164-3         160-6 

do 

Determine  as  carefully  as  j-ou  can  the  value  of  ^  when  C=2-6  hours. 

-MISCELLANEOUS  EXAMPLES. 

1.  Find  a  polynomial  of  the  third  degree  in  x  which  assumes  the  values 
37,  11,  13,  29  when  x  is  equal  to  2,  1,  -1.  -2  respectively.  Also  find  the 
minimum  value  of  the  function. 

2.  A  steam  electric  generator  on  three  long  trials,  each  with  a  different 
point  of  cut  ofl"  on  steady  load,  is  found  to  use  the  following  amounts  of  steam 
per  hour  for  the  following  amounts  of  power  :  — 

Lb.  of  steam  per  hour  40-20  6650  10800 

Indicated  horse-power  210  480  706 

Kilowatts  produced  114  290  435 

Find  the  indicated  horsepower  and  the  weight  of  steam  used  per  hour  when 

330  kilowatts  are  being  produced. 

3.  Using  first  differences  only,  find  the  value  of  ./■(4817-36)  if 

,r  4S16  4817  4818 

./•(.(■)         693974063         694046108         694118146 
To  what  extent  may  the  result  be  relied  on  ? 

4.  Interpolate  for  .i-  =  35"2967  from  the  following  table  and  state  the 
accuracy  of  the  result  :- 


.,• 

/(•'•) 

35-28 

I6202I -40139 

35-29 

162043-22198 

.35-30 

162065  04777 

5.  If  the  sun's  apparent  longitude  at  O.M.N,  on  1914  Aug.  1st  was 
128°  22'  25"  8  ;  on  Aug.  3rd  130°  17'  15"  2  ;  on  Aug.  5th  13-2'  12'  7  "8  ;  find 
the  longitude  at  G.M.N,  on  Aug.  2nd. 

6.  At  the  following  draughts  in  sea  water  a  particular  vessel  has  the 
following  displacements  : — 

Draught  in  feet  -         -  15  12  9  6 

Displacemeat  in  tons  -         "2098         1512  1018  536 

What  are  the  probable  displacements  wheo  the  draughts  are  11  and   13  feet 

respectively  ? 

7.  The  following  table  gives  the  e.xpejtatioa  of  life  for  males  of  the  ages 
mentioned  : — 

Present  age         ....         5  10  15  25  35 

Expectation  of  life  in  years       -      5087      47  60       43-41       35-68       28  64 

Find  the  expectations  of  3  men  whose  ages  are  18,  20,  30  years  respectively. 

o.  3 
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8.  The  H*'  table  of  the  Institute  of  Actuaries  gives  the  following  values 
of  h  :- 


a- 

^ 

10 

100000 

12 

99113 

14 

98496 

16 

97942 

18 

97245 

20 

96223 

Obtain  the  values  of  Ix  when  x  has  the  values  11,  13,  15,  17,  19. 

9.  The  sun's  apparent  right  ascension  at  mean  noon  is  given  in  the  follow- 
ing table  : — 


,  June  \ 

4 

34 

0-69 

6 

4 

54 

31-75 

11 

5 

15 

10-64 

16 

5 

35 

55-34 

21 

5 

56 

43-18 

26 

6 

17 

30-99 

Determine  the  apparent  right  ascension  for  1914  June  15d.  16h. 

10.  Let  the  number  of  recruits  in  the  United  States  army  whose  height 
did  not  exceed  x  inches  be  denoted  by  y.  Determine  the  number  whose 
height  exceeded  66  inches,  but  did  not  exceed  66J  inches,  if 


62 

618 

63 

1855 

64 

3802 

65 

6821 

66 

10296 

67 

14350 

68 

17981 

69 

21114 

70 

23189 

71 

24674 

sun's  apparent  declination  3 

at  mean  noon  i 

1914,  May  1 

+  14° 

54'    19" -6 

3 

+  15° 

30'    25" -0 

5 

+  16° 

5'    28" -4 

7 

+  16° 

39'    27" -6 

9 

+  17° 

12'    20" -3 

11 

+  17° 

44'    4" -3 

Determine  the  hourly  difiference  in  9  for  1914  May  5  and  7. 


i 


0-3 

0-4 

0-5 

0-6 

0-7 

O'S 

0-9 

1-0 

1-1 

0-2S1 

0-525 

0-S54 

1-2G0 

1-725 

2-2-29 

2-753 

3-2S2 

3-Sll 
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12.  In  the  following  table  •■*  denotes  the  space  described  at  the  time  t  by 
a  point  in  a  mechanism.  Tabulate  the  values  of  the  velocity  and  the  acceler- 
ation, and  obtain  their  values  when  <  =  0-45. 

t  (time  in  sees.)       0      0-1 
I  (space  in  feet)       0    0-03-' 

13.  There  is  a  piece  of  mechanism  whose  weight  is  200  lbs.  The  following 
values  of  s  in  feet  show  the  distance  of  its  centre  of  gravity  from  some  point 
in  its  straight  path  at  the  time  t  seconds  from  some  era  of  reckoning.  Find 
its  acceleration  at  the  time  ^=2-05,  and  the  force  in  pounds  which  is  giving 
this  acceleration  to  it. 

t       2-00  2-02  2-04  2-06  2 -OS  2-10 

s     0-3090        0-4931        0-6799        0-8701         10643         1-2631 

14.  The  following  values  of  p  (lbs.  per  sq.  foot)  and  ^^C  are  for  saturated 
steam.     Calculate  with  as  much  accuracy  as  the  numbers  will  allow  the  value 


dp 
of  ^fore=115'C. 

d 

105 

110 

115 

1-20 

125 

V 

2524 

2994 

3534 

4152 

4854 

15.   From  the  following  table  of  the  secular  variations  of  the  magnetic 
declination  determine  the  chief  period  : — 

Year.  Declination. 

1540-5  -    7" -463 

1580-5  -   8° -500 

1620-5  -   6° -250 

1660-5  -   0°-500 

1700-5  +   8=  "2-28 

1740-5  +15°-755 

1780-5  +20-832 

18-20-5  +2-2° -380 

1860-5  +19° -364 

19O0-5  +14^-833 
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CHAPTER    III 

THE   CONSTRUCTION    AND    USE   OF 
MATHEMATICAL   TABLES 

13.   Interpolation  by  First  Differences. 

All  mathematical  tables  are  alike  in  that  they  give  numerical 
values  of  the  function  for  certain  values  of  the  argument,  which 
are  so  chosen  that  intermediate  values  of  the  function  may 
be  derived  by  interpolation.  The  interval  of  the  argument 
varies  in  the  different  mathematical  tables,  its  choice  in 
each  case  being  determined  by  the  rapidity  of  variation  of  the 
function.  In  most  of  the  published  elementary  tables,  vi^hich  are 
intended  for  the  use  of  students,  the  interval  of  the  argument  is 
so  chosen  that  the  second  difference  is  smaller  than  one  unit  in 
the  last  place  of  decimals  retained :  on  this  account  it  is  not 
necessary  to  consider  second  differences  in  the  interpolation,  and 
the  function  can  be  calculated  for  any  value  of  the  argument 
intermediate  between  two  of  the  tabulated  values  by  a  simple 
proportion  in  the  following  manner  : — 

Consider  a  function  fix)  which  is  tabulated  for  the  values 
a*o,  Xq  +  w,  Xq  +  2w,  ...  of  the  argument.  Then,  since  the  function 
increases  uniformly,  when  digits  beyond  a  certain  place  of  decimals 
are  neglected,  we  have 

/(■tq  +  w')  -/(xq)  _  vy 

/{Xo  +  W)-/{Xo)    "     w 

where  0<tv'<w. 

From  this  it  follows  that 

w' 
/(a^o  +  W)  =/(a-o)  +  —  {/{Xq  +  w)  -/{xq)  }, 

so  that,  if  the  first  difference  /{xq  +  w) -/(x^)  is  known,  the  value 
of  the  function  for  any  value  of  x  between  x^  and  a-^  + «"  may  be 
found.  This  equation  is  known  as  the  "  Rule  of  Proportional 
Parts."  Geometricalh^,  what  we  do  is  to  replace  the  curve  through 
two  consecutive  points  by  the  chord  joining  them. 
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If  the  first  differences  are  only  sensibly  constant,  the  result 
obtained  will  be  subject  to  a  certain  error.  If  w'  =  n  w  where 
0<n<  1,  then  the  error  may  be  put  in  the  approximate  form 

n(n-l)  , 
tu-J    {x  +  nw) 

where  0<?i'<  1. 

Example  1. — In   the    following  table    the  entry  is  the  cube  root  of  the 
corresponding  argument.     Find  the  cube  root  of  1619"25. 
Argument.  Entry.  A 

1619  11 •7421858 

24171 

1620  11-7446029 

24161 

1621  11-7470190 

If  we  desire  accuracy  up  to  the  fifth  decimal  place,  then  we  may  regard  the 
first  differences  as  constant. 

Here  .fo=1619     tr=l     if'  =  0-25. 

Hence      v'lG19"25=  11-74219  +  jnfo  (0-00242) 
=  11-74280. 

Example  2.  —  If  a- =  1000  show  that,  when  second  differences  are 
neglected,  the  error  committed  in  computing  by  this  method  the  value  of 
log,(,(1000  +  "-),  where  0  <  ic  <=  1,  will  have  no  effect  on  the  sixth  decimal 
place  of   the  logarithm. 

Here  lug  (x  +  v)  =  log  x  +  v  { log  (./;  -i- 1 )  -  lug  x  \ 

and  the  error  committed  is 

E="i^^r{x-,vy) 

where  0 -=  w' -=  1     and    /(.i')  =  log]o,f. 

I  I        I  f^"-'  I 

\J"{X+W')       =1^2    logioU'  +  w') 


Now 


I  {x  +  w'f  I 
where     m  =  modulus  of  common  logarithms  =  0  43429. 
But  the  maximum  value  of  io(\  -  )/■)  is  \  since 

w(l-(6-)  =  i   -(^-,,f  <i. 

Hence  the  error  committed  is  lees  than 

0-43429 


E*c 


(1000 +  ('.-')- 

1 


16 .  10'' 

0  0000000625, 
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which  shows  that  the  error  committed  does  not  affect  the  sixth  decimal  place 
in  the  logarithm  of  the  number  (1000+  w). 

Example  ^1— The  differences  between  the  six-tigure  logarithms,  to  base  10, 
of  the  successive  integers  between  5000  and  5100  are  sensibly  constant. 
Prove  this  and  determine  the  value  of  the  constant  difference. 


14.    Case  of  Irregular  Differences. 

When  the  second  difierences  cannot  be  neglected,  however, 
this  method  cannot  be  applied.  This  happens  frequently  when 
the  rate  of  change  of  the  function  is  large  compared  with  that  of 
the  argument,  as  is  the  case  with  the  logarithmic  sines  and 
tangents  of  small  angles.  We  must  then  use  the  formulae  of 
Chapter  II.,  including  differences  to  the  order  required. 

Example  1. — To  obtain  log  sin  0°  16'  24" "5. 
From  Schron's  tables  we  have 

A  A- 

i/ sin  0°  16' 10"  =  7  ■6723450 


A^' 


iy  sin  0°  16' 20"  =  7 -6767993  -452 

44091 


i/ sin  OM  6' 30"  =  7 -68 1208  4 

-443 

i>sinO°16'40"  =  7-6855732 

43648 

Using  Gauss'  formula 

n{n 
f(a  +  n  w)  - /;,  +  n5J^+      ^ 

-1)^„,       {n  +  l)n{>i-l) 
!        ^'-^o  +             3  ! 

obtain 

iysin0°16'24" 

•5-7-6767993-1 

19841 

56 

=  7-6787889 

This  is  the  answer  correct  to  7  places,  the  value  correct  to  10  places  being 
found  from  Vega  to  be  7-6787889383. 

The  value  obtained  by  the  "  Rule  of  Proportional  Parts"  is  7-6787834. 
This  result  is  correct  only  to  4  places  of  decimals,  and  shows  the  necessity  of 
interpolating  by  some  such  method  as  the  above  the  value  of  the  logarithmic 
sine  of  a  small  angle.  This  may  be  shown  to  be  also  true  in  the  case  of  the 
logarithmic  tangents  of  small  angles. 

Example  ;^.  — Compute  the  value  of  log  tan  0°  15'  35". 
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15.  Construction  of  New  Tables. 

In  spite  of  the  large  number  of  different  mathematical  tables 
that  have  been  constructed  to  assist  the  computer,  further 
research  in  the  various  branches  of  mathematics  and  the  kindred 
sciences  shows  that  these  are  not  sufficient  for  all  purposes.  A 
research  student  often  requires  new  tables  for  the  question  which 
he  is  discussing,  and  these  he  must  construct  for  himself. 

The  first  question  to  be  settled  is  the  degree  of  accuracy  to 
which  the  tables  are  to  be  computed,  i.e.,  the  number  of  decimal 
places  to  be  retained.  This  varies  in  different  cases ;  but  it  is 
always  desirable  to  extend  the  original  calculations  beyond  the 
number  of  places  which  are  ultimately  to  be  retained,  in  order  to 
avoid  the  cumulative  effect  of  the  neglected  digits. 

The  general  plan  of  the  calculation  is  to  compute  first  with  great 
accuracy  the  values  of  the  function  for  a  series  of  widely-spaced 
values  of  the  argument.  For  this  purpose  infinite  series  are  generally 
employed.  The  values  of  the  function  for  a  more  closely-spaced 
series  of  values  of  the  argument  are  then  derived  from  these  by 
means  of  the  formulae  of  Chapter  II.  :  this  second  process  is  called 
Suhtabulation.     We  shall  consider  these  processes  in  §§  17,  18. 

16.  Tabulation  of  a  Polynomial. 

The  simplest  case  to  be  considered  is  the  tabulation  of  a  polynomial. 

If  the  degree  of  f{x)  is  n,  it  will  not  be  necessary  to  compute  its  value 
directly  for  more  than  (n-l-1)  values  of  the  variable,  as  may  be  seen  from  the 
following  examples. 

Example  1.  —Let    f{x) -x^  +  Zx--\-2x-\. 

When  .1=0,  1,  2,  3,  we  have  /(.t)=  -1.  5.  23,  59.  We  thus  obtain  the 
following  table  :  — 

X  J\x)  A  A-  A' 


3G 

3  5<J 

Since  J'(x)  is  of  the  third  degree  in  x,  its  third  differences  will  be  constant ; 

in  this  case  they  will  be  equal  to  G.     Extending  tiie  column  A'  by  inserting 

6'8  we  can  then  extend  the  column  A-'  by  simple  additions  or  subtractions. 

A  repetition  of  the  same  process  will  enable  us  to  extend  the  column  A,  and 
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thereafter  to  determine  the  values  of  /(.r)  for  all  integral  values  of  x.     The 
resulting  table  is 


X 

./■{■<■] 

A 

A- 

-3 

-7 

6 

-12 

-2 

-1 

0 

-6 

-1 

-1 

0 

0 

0 

-1 

6 

6 

1 

5 

18 

12 

2 

23 

IS 

3  59  24 

60  6 

4  119  30 

90  6 

5  209  36 

By  this  simple  method  we  are  enabled  to  determine  the  values  of  a  poly- 
nomial fix)  for  all  positive  and  negative  integral  values  of  x.  We  can  then 
use  any  of  the  formulae  of  interpolation  already  established  to  obtain  the 
values  of  /(x)  for  fractional  values  of  x.  The  method  may  be  extended  to 
functions  which  are  not  rational  or  integral,  but  in  this  case  it  will  be 
necessary  to  make  the  interval  between  successive  values  of  x  so  small  that 
the  differences  of  some  definite  order  are  constant  throughout  the  range  of 
values  under  consideration. 

Example  3.— Form,  the  table  of  values  of  /(.r)^.r'*  +  5.f^- 7.i"-f.r-- 4  for  the 
values  0,  1,  2,  3,  4  of  the  argument  x,  and  then  deduce  the  values  of  J\x) 
corresponding  to  the  values  -  1 ,   -  2,   -  3,  -  4  of  x. 

17-  Calculation  of  the  Fundamental  Values  for  a 
Table  of  Logarithms. 

We  shall  now  illustrate  the  process  of  computing  tables  of 
functions  other  than  polynomials  by  describing  the  calculation 
of  a  table  of  logarithms. 

The  first  stage  in  this  consists  in  determining  the  logarithms 
of  all  the  prime  numbers  between  certain  limits.  It  is  evident 
that  when  the  logarithms  of  the  primes  are  known,  the  logarithms 
of  all  composite  numbers  can  be  derived  from  them  by  mere  use  of 
the  formula 

log  a  6  =  log  a  +  log  b. 
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In  order  to  obtain  the  logarithms  of  the  primes  we  take  two 
equations  whose  roots  are  integers  and  which  differ  only  in  their 
constant  term.      Such  a  pair  are 

or*  -  3  X  -  2  =  0. 

Now  x'-3x  +  2  =  {x-\)- {x  +  2) 

oc'-3x-2  =  {x+l)-  {x-2). 
Hence 

lug  (x  -\)  +  h  log  {X  +  2)  -  log  (a;  +  1 )  -  1  log  (.;  -  2) 

,  .      o;^'  -  3  a-  +  2 


a  '"o 

X"  - 

1  + 

3 

x-2 
2 

Hog 

x^ 

-3x 

1  - 

2 

3x 

which,  by  use  of  the  ordinary  logaritiunic  scries,  may  be  written, 
if  the  logarithms  are  all  to  the  base  e, 

x'-3x       ''  Xx'-^x)        '  yx'-3xf       '" 
This  is  called  Borda's  Formula. 

Putting  x  =  b,  6,  7,  8,  we  obtain  the  four  equations 

log  2  -  \  log  3  +  h  log  7  =  0-0181838220854376... 

ilog2  +    logo-     log7=0-0101013536587597... 

-  2  log  2  +  2  log  3-\  log  5  =0  0062 11 2599992784... 

-  e  log  3  +  ^  log  5  +    log  7  =  0-0040983836020895... 

Solving  these  four  equations  we  obtain  the  logarithms  (to  the 

base  e)  of  the  prime  numbers  2,  3,  5,  7,  viz.  : — 

log  2  =  0-6931471805599453... 
log  3  =  1-09S6122886G81096... 
log  5  =  1-6094379124341003... 
log  7  =  1-9459101490553133... 

Further  substitution  for  ./;  in  Borda's  formula  will  give  immetli- 
ately  the  logarithms  of  the  other  prime  numbers.  For  example,  on 
putting  X  =  9  the  left-hand  side  becomes 

log8+Uogll-loglO-'  l„g7. 


42  INTERPOLATION  AND  NUMERICAL  INTEGRATION     [ch.  hi 

Of  these  log  7  has  been  found ; 

10^8    -31og2  =  2-0794415416798359...  ; 
and  log  10  =  log  2  +  log  5  ^  2-3025850929940456.... 

Hence  we  obtain  log  11  =  2-3978952727983705  .... 

It  is  to  be  noted  that  these  logarithms  are  to  the  base  e.  But 
as  we  have  now  found  log,  10,  we  can  at  once  derive  its  reciprocal 
logjoe  =  0-4342944819032518...,  and  all  the  above  logarithms 
can  now  be  converted  to  the  base  10  by  multiplying  them  by  this 
"  modulus." 

We  thus  obtain  the  table 

logio  2  =  0-3010299956639812 

logio  3  =  0-4771212547196624 

logio  5  =  0-6989700043360188 

logjo  7=0-8450980400142568 

logroll  =  1-0413926851582250 

etc. 

For  the  calculation  of  the  logarithms  of  the  larger  primes  Haros'  Formula 

will  be  found  more  useful.     In  this  case  consider  the  equations 

a;*-25a;2+144  =  0. 
From  these  we  deduce 

{x  +  5)  (.r  -  5)  x"- xi-25x"-       _         (.c^  -  25  x-^- +  72)  -  72 

^°g  (x  +  ^)(x-Z)[x  +  ^)\x  -  4)  -  ^°S  x'  -  25 a;2  +  144  "  '°§  [x^  -  25 x^-  +  72)  +  72 

whence,  if  the  logarithms  are  to  the  base  10,  we  have 

log{;r  +  5)  =  {log(.c  +  3)  +  log(,f-3)  +  log(.r  +  4)  +  log(.r-4)}-{log(.r-5)+21og.r} 

--"^  {.^^251^2  +  *  (.^25'a?T72)  +  "  / 
where  m  =  logif,  e. 

Example  Z.— To  obtain  the  logarithm  of  the  prime  number  43  to  the 
base  10. 

Let  .1  +  5  =  43. 

Then  >i  =  38 

log  43  =  log  4 1  +  log  35  +  log  42  +  log  34  -  log  33  -  2  log  38 
-  0-8685889638065036  {0-0000351372402040  +  0-0000000000000145}. 
But  the  logarithms  of  41,  35,  42,  34,  33,  38  are  assumed  to  have  been 
previously  calculated.  Hence 

log  43=  1  -6127838567197355  -  1  -5185139398778875 
1  -544068044350-2756  3-1595671932336203 
I  -6232492903979005  0  0000305 1 98 1 90724 
1-531478917042-2551 
=  1-6334684555795865 
Example  2.— Show  that  logio  131  =  2-1172712956557643. 
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18.    Amplification  of  the  Table  by  Subtabulation. 

Having  now  obtained  the  logarithuKs  of  the  prime  numbers 
within  a  certain  range,  we  derive  from  them,  by  mere  addition,  the 
logarithms  of  the  composite  numbers  which  lie  in  this  range,  and 
thus  obtain  the  logarithms  of  all  the  integers  within  the  range. 
We  shall  now  show  how,  by  subtabulation,  the  extent  of  this  table 
may  be  increased  manyfold.  It  will  be  supposed  that  7-place 
accuracy  is  required. 

For  example,  suppose  the  logarithms  of  the  integers  31,<32,  33, 
34,  35,  36  have  been  obtained  directly  ;  we  shall  show  how  the 
logarithms  of  all  the  integers  between  330  and  340  can  be  derived 
from  them. 

Neglecting    the    characteristics   of    the    logarithms,   which   can 
always    be   written    down    by   inspection,    we  form    the   following 
table  : — 
X  logx-  A  A-  A3  A* 

31  -49136169 

137^829 

32  -50514998  -42433 

1336396  2535 

33  -51851394  -39898  -223 

1296498  2312   (-208) 

34  -53147892  -37586  -192 

1258912  2120 

35  -54406804  -  35466 

1223446 

36  -55630250 

Putting  a  =  33,   h  =  ^l,   iv  =  1  in  Gauss'  formula 

ri(l-?t),,  .       n(l-n')       . 
J(a  +  n  IV)  =Jo  +  'i  oJ,^ Y\ —  ^'-^^ 3~l *^i 

^n(l-.0(2-rO       , 
4  !  '     -^  i 

we  obtain 

y'(33-l)=7;  +  0-l  0/'  -0-045  oVo  -U-Ul65o7J  + 0-0081  iio'j\^ 
=  -51851394-38 
129650       1 
1795 
=  -51982800. 
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Thus  to  7  decimal  places  log  331  =  2-5198280.  But  it  is  not 
necessary  to  calculate  the  values  for  all  the  numbers  between 
330  and  340  in  this  way.  It  is,  in  fact,  simpler  to  calculate 
log  332  from  the  formula 

f{a  +  -^  to)  =/(«  +  -Ay  iv)  +  {/{a  +  j-V  tv)  -/(a) }  +  8, 
where     &,  =/(«  +  f^  tv)  -  2/ (a  +  j\  w)  +/ (a) 

=  0-01  syo  +  0-001  sy^  -  o-ooi  /.  sy, 

and  then  to  calculate  log  333  from  the  formula 

/(«  +  TO  ^)  =/(«  +  TO  ^)  +  {/(«  +  TO  «0  -/(«  +  TO  W') }  +  Sa 

where     S^  =f{a  +  ^^  w)  ~  2/ (a  +  -f-^  to)  +f{a  +  ^\  w) 

=  001  sy  +  0-002  sy,  -  00017  /t  sy; . 

The  complete  table  of  these  quantities  S  is  as  follows  : — 

k  =  0-01 8y  +  0-001 8y,  -  o-ooio  /x  sy; 
^3  =  0-01  sy;  +  0-002  sy  -  0-0017  /.  sy; 
S4  =  0-01  s^/;  +  0-003  8^/;  -  0-0019 IX  sy; 
Sg  =  0-01 3y;  +  0-004  sy;  -  0-0021  /^  sy; 
Sg  =  0-01  sy;  +  0-005  sv;  -  0-0020  ix  sy 
s,  =  0-01  sy  + 0-006  sy; -00021  /xsy; 
S3  =  0-01  sy  +  0-007  sy;  -  0-0018  /.  sy; 
Sg  =  0-01  sy;  +  0-008  sy;  -  0-0017  /^  sy; 
Sio  =  0-01  sy;  +  0-009  sy-  -  0-001 2  /.  sy- 

The  numerical  work  is  arranged  as  in  the  following  table,  in  which 
the  logarithms  required  are  found  in  the  last  column. 


-  ; 

39898-0  +  2 

312-0-2 

108-0 

129649-8 

2-51851394 

To 

+ 

1795-4- 

38-1- 

1-7=  - 

-1755-6 

131405-4 

2-519827994 

TO 

- 

399-0  + 

2-3  + 

0-2=  - 

-  396-5 

131008-9 

2-521138083 

:3 

Tcr 

- 

399-0  + 

4'6  + 

0-4=  - 

-  394-0 

130614-9 

2-522444232 

4 
To 

- 

399-0  + 

6-9  + 

0-4=  - 

-  391-7 

130223-2 

2523746464 

To 

- 

399-0  + 

9-2  + 

0-4=  - 

-  389-4 

129833-8 

2-525044802 

6 

To 

- 

399-0  + 

11-6  + 

0-4=  - 

-  387-0 

129446-8 

2-526339270 

To 

- 

399-0  + 

13-9  + 

0-4=  - 

-  384-7 

129062-1 

2-527629891 

8 

To 

- 

399-0  + 

16-2  + 

0-4=  - 

-  382-4 

128679-7 

2-528916688 

i. 

- 

399-0  + 

18-5  + 

0-4=  - 

-  380-1 

128299-6 

2-530199684 

10 
T^ 

- 

399-0  + 

20-8  + 

0'2=  - 

-  378-0 

127921-6 

2-531478900 
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The  method  followed  in  constructing  this  table  will  best  be  shown  by 
considering  two  examples. 

In  the  line   I'u 

+  1795  4  is  the  numerical  value  obtained  for     0'045    S^/o 

-  38-1  ,,  ,,  „  -0-0165  Sy^ 

-  1-7  ,,  >,  ,,  O'OOS    fi 5*y,^ 
1755 '6  is  the  sum  of  the  three  preceding  numbers. 

129649-8  is  the  numerical  value  obtained  for  001  5/^  . 
131405-4  is  the  sum   of   the   two   preceding   numbers,   and   therefore 
represents 

001  6/,  +  0-045  5-/u  -  0  0165  Pj\  +  0-008  m  5\r, . 
2-51851394  is  the  numerical  value  of  /„ 

2-519327994,  which  is  obtained  as   the  sum  of  the  two  preceding 
numbers,  represents 

/a  +  0-01  5/j^  +  0  045  oVo  -  0-0165  5^J\  +0-008  n  sy^  • 
In  the  line  ^^ 

-  399-0  is  the  numerical  value  obtained  for     0-01    5-/o 
+     2-3  „  „  „  0-001 6Y, 

+     0-2  ,,  ,,  ,,  -0-001  A*  5V', 

-.396-5,  which  is  obtained  by  summing  the  three  preceding  numbers, 
represents 

0-01 52/0 + 0-001  av,  -  0-001  ^  sy^  or  s.,. 

131405-4,  which  has  previously  been  obtained,  represents 

1310US-9,   which  is  obtained  by  summing  the  numbers    immediately 
above  it  and  to  the  left  of  it,  represents 

Now  2-519827794,  as  we  have  already  seen,  represents  f(a  +  -^w). 

2521138083,  which  is  obtained  by  adding  the  numbers  immediately 
above  it  and  to  the  left  of  it,  represents 

/(a+.'ir"-). 

The  results  in  the  last  column  are  the  logarithms  of  the  numbers 
330,  331,  ...  340.  They  cannot,  of  course,  be  relied  on  up  to  the 
last  figure,  but  will  be  found  correct  up  to  seven  decimal  places. 
This  accuracy  may  be  obtained  by  taking  one  decimal  more  in  the 
first  difference  than  in  the  tabular  function,  one  decimal  more  in 
the  second  difference  than  in  the  first,  and  so  on,  when  these  can 
be  obtained.  A  check  on  the  accuracy  is  obtained  by  carrying 
on  the  process  until  log  340  is  computed.  Comparison  of  this 
result  with  the  value  for  log  34  will  show  whether  the  computa- 
tion has  been  accurately  performed. 
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We  may  remark  that  it  is  an  improvement  on  the  foregoing 
process  if  the  differences  required  for  the  subtabulation  are 
computed  independently  from  the  expressions  in  series,  instead  of 
being  derived  by  forming  a  table  of  differences. 

Example  i.— Given  log  24  =  1  -3802112 

log  25  =  1-3979400 
log  26  =  1-4149733 
log  27  =  1-4313638 
log  28  =  1-4471580 
log  29  =  1-4623980 
obtain  the  logarithms  of  the  numbers  261,  262,  ...  269. 

19.    Radix  Method  of  Calculating  Logarithms. 

When  the  logarithm  of  a  number  is  required  to  a  large  number 
of  decimal  places,  its  determination  either  by  means  of  an  infinite 
series  or  by  interpolation  becomes  very  laborious  and  liable  to 
error.  A  much  simpler  and  more  useful  method  in  such  a  case  is 
that  known  as  the  Radix  Method. 

Let  N  be  the  given  number  whose  logarithm  is  required  to  n 
places  of  decimals.     Then  since 

#=10^.a.(l+iV^o) 
where  N^^  is  a  decimal,  x  an  integer,  and  a  =  0,  1,  2,  ...  9, 

logio  N=x-v  log,o a  +  logio (1  +  N^. 
Thus,  to  determine  logioiV,  it  will  be  sufficient  to  compute  logma 
and  log,o(l  +  ^o)-    The  foniier  has  been  obtained  above.     We  shall 
now  show  how  to  evaluate  the  lattei-. 

Let  rj  denote  the  decimal  consisting  of  the  first  m  figures  in 
iVo,  and  let  N^  be  the  difference  of  N^  and  r^.     Then 

(l  +  iV^„)-(l+rO  =  i\^P 

Now  divide  N^  by  (l+?'i)  until  there  are  m  figures  in  the 
quotient  r^.     If  we  call  the  remainder  iV^,  then 

Similarly,  if  the  remainder,  after  dividing  iV^o  by  (l+ri)(l  +  r2) 
until  there  are  m  figures  in  the  quotient  j-,,  is  denoted  by  N^^  we 
have 

and  in  general, 

N^_,  -  (1  +r,)  (1+r,)  ...  (1  +r^_i)r^  =  iV;. 
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These  equations  may  be  written  in  the  form 

(1+iVo)  -.V,    =(i+n) 

iV^j      -i\r,     =(l+r,)(l+r,)-(l  +  r,) 

JVo     -N.     ^{l+r,)(\+r^{l+r.^-{\+r,){l+r„) 

iV^,-.-^p-i  =  (l+n)(l+^--..)---(l+'V-i) 

-(l+n)(]+r,)...(l+7v_.) 
^p-,-^p     =(l+r,)(l+r,)...(l+0 

-(l+rO(l+n)...(l+r,_,). 
Adding,  we  obtain 

(l+JYo)    -N^     =(l+rO(l+r,)...(l+r„) 

whence 

1  +.^-0     ={l+r,)  (l+r,)  ...  (1  +r,)  +  xAV 

If,  therefore,  the  remainder  A^^,  has  no   significant   figures    in    its 
first  n  decimal  places, 

log{l+ayo)    =log(l+r,)  +  log(l+r,)+  ...  +log(l+r,J. 
The  quantities  (l+rj),  (l+r.,),  ...  are  called  the  radices.     As 
they  will  be  of  the  form  (l+  --^  j  where  k  and  I  are  integers,  and 

as  k  is  always  small  compared  with  10'  we  may  compute  log  (1  +  r^) 
by  means  of  the  logarithmic  expansion 

log  (1  +  7^)  =  r^-h  r/  +  i  r^'  -  ^  r/  +  . . . 
The  values  of  the  logarithms  of  these  radices  when  k  =  0,  1,  2,  ... 
999,  and  Z  =  3,  6,  9,  12,  15,  18,  21,  24,  as  well  as  of  the  numbers 
1,  2,  ...  9,  have  been  computed  to  twenty-four  places  of  decimals.* 

The  inverse  process  of  determining  the  antilogarithm  corre- 
sponding to  a  given  logarithm  is  equally  simple.  When  the 
logarithm  is  not  in  the  table,  we  have  merely  to  take  out  the 
next  lower  and  subtract  it  from  the  given  one.  Then  take  out 
the  next  lower  and  subtract  it  from  this  remainder ;  and  so  on 
until  the  remainder  consists  of  zeros.  The  product  of  the  radices 
corresponding  to  the  logarithms  taken  out  is  the  number  required. 

*  Gray,  Tables  for  the  formation  of  Logarithms  and  Antilogarithms  to 
twenty-four  or  any  less  number  of  places.     London,  C  <fe  E.  Layton  (1900). 
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We    sliall     illustrate    these    two    methods    in    the     following 
examples : — 

Example  /.—To  determine  logjo43  to  24  places  of  decimals. 

Here  .^=1,  a  =  3. 

Taking  to  =  3,  as  in  Gray's  tables,  we  have 

1  +  xVo  =  1  •  433  333  333  333  333  333  333  333 
l  +  ri  =  l-433 
iVj  =  0  •  000  333  333  333  333  333  333  333 
r.3=0-000232 

iV2=0'  000  000  877  333  333  333 .333  333 
(l+rj)/-2  =  0- 000332456     i.e.,  0-000  333.333  -  0-000000  877 
(I+ri)(l  +  r2^  =  (l-fr,)  +  (l+ri)r, 
=  1-433.332  456. 
The  above  will  be  sufficient  to  indicate  the  general  theory.     The  process 
can,  however,  be  simplified  by  arranging  the  work  in  the  following  manner  :— 
1  -  433  333  333  333  333  333  333  333 
1-433 


Diff.=  0-000  333  333  333  333  333  333  333 

1  -  433  )  0-  000 333 .333 ,333  333  333  333 .3.33  (  0-  000 2.32 
286  6 

46  73 
42  99 

3  743 

2  866 

877  =  Remainder. 
1-4.33333  333 

0-000000  877 
Diff.  =  1-433  332  456  =  New  Divisor. 

1  •  433  332 456  )  0-  000 000  877  333 333 333  333  333  (  0-  (000)2612 
859  9994736 

17  333  859  73 
14  3333-24  56 


3  000535173 
2  866  664  912 


133  870261  =  Remainder. 
1  -  433  333  333  333  333  333 
0-000  000000133  870261 
Diff.  =  1-4.33333  333  199  463  072"=  New  Divisor. 
Should  the  original  dividend  not  be  known   beyond   the   24th   place   of 
decimals,  we  should  now  have  to  proceed  by  a  method  of  contracted  division. 
In  the  present  case  we  know  the  figures  after  the  24th  decimal,  and  may 
therefore  continue  as  above.     But  we  shall  adopt  the  contracted  method. 

*  0- (000)2  612  is  written  in  place  of  0-0000C0612 
0(000)3  093  „  „  0-000000000093  and  so  on. 
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The  new  divisor  contains  19  figures  ;  the  new  dividend 
■(000)»  133  870  261  383  833 
contains  only  15  significant  figures  after  the  point.     We  must  therefore  cut 
off  4  figures  from  the  divisor.     Thus  we  obtain 

1  •  433  333  333  199,4,6  )  0"  (000)^  133  870  -^61  333  333  (  0"  (000 )3 093  * 
128  999  999  987  951 


4  870  261345  382 
4  299  999  999  598 


570261  345  784  ::^  Remainder. 
1-433  333  333  333 333  333  333  333 

0-  000  000  OUO  000  570  261  345  784 


DifiF.  =  1-433  333  333  332  763071  987  549  =  New  Divisor. 
As  the  new  dividend  0"  (000)^  570  261  345  784  contains  only  12  significant 
figures  and  the  divisor  -25  we  must  cut  oti'  13  figures  from  the  latter.  It 
then  becomes  1-433  333  333  33.  As  these  are  the  figures  in  the  original 
dividend  there  will  be  no  further  new  divisors.  This  will  usually  happen 
when  the  fourth  division  is  reached.  We  may  therefore  determine  the 
remaining  figures  by  continued  contracted  division. 

1  -  433 333  333  33  )  0'  (000)^  570 261  345  784  (  0 •  (000)»  397  856  752  873 
430000000  000 


140-261345  784 
129  000000000 

11261345  784 
10033  333  333 

12-28  012451 
1146  666  666 

81345  785 
71666  667 


9  679]  18 
8  600  000 

1079118 
1  003  333 


75  785 
71667 
Tll8 
2867 

1251 
1146 

105 
100 

5 
4 
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The  remaining  radices  are  thus 
1  ■  (000)^  397 
1- (000)^856 
l-(000)«752 
1- (000)7  873 

Using  Gray's  tables  of  the  radices  we  have 


log  3 

0-  477  121  254  719  662  437  295  028 

log  1-433 

0- 156  246  190  397  344  475  994  693 

log  1- 000232   = 

100  744  633875  845  680  796 

log  1- (000)- 612  = 

265  788141593627  087 

log  1- (000)3  093 - 

40  389  386  815124 

log  1- (000)^  397 - 

172414  909316 

log  1- (000)5  856  = 

371756077 

log  1- (000)6  752  = 

326589 

log  1- (000)'' 873  = 

379 

Sum  = 

:  0-633  468  455  579  586  526  405  089 

Hence  since  x-l  the  required  1 

ogarithm  is 

log  43  =  1  -  633  468  455  579  586  526  405  089. 

Example  2. — To  determine  the  antilogarithm  of 

1  •  120  619  882  724  133  396  646  450. 
As  the  integral  part  only  affects  the  characteristic  we  shall  neglect  it  for 
the  present. 

Using  Gray's  tables  we  then  obtain 

0- 120  619  882  724  133  396  646  450 
log    1-320  =0-120  573  931205  849  868  472  706 


Diff. 

45951518  283  528173  744 

log  1-000105   = 

45  598  5-26  719080137349 

Diff. 

352  991564448  036  395 

log  1 -(000)^  812  = 

352  646  976130  787  551 

Diff. 

344  588  317  248  844 

log  l-(000)-^793  = 

344  395  524  012  7-26 

Difl'. 

192  793-236118 

log  1  (000)^43  = 

192  392  455  483 

Diff. 

400  780635 

log  1 -(000)-^  922  = 

400419512 

Diflf. 

361  123 

log  1-(000)«831  = 

360899 

Diff. 

224 

log  1 -(000)' 516  = 

224 

Diff. 

0 

Multiplying  these  radices  together  we  have 
1 -(000)^43  X  1 -(000)3  922  X  l-(000)«  831  ^  1.(000)7  516. 


1  (000)^  443  922  831  516. 


I 


19.  20]  MATHEMATICAL  TABLES  51 

This  product  of   the  last  four  factors  can  always  be  written  down  by 
inspection.     The  other  multiplications  may  be  arranged  as  follows  : — 
1  •  000  000  000  000  443  922  831  516 
1-000000000  793 


1  •  000  000  000  000  443  922  831  516 

700  000000  000  310 

90000000  000040 

3  000000  000000 

1  •  000  000  000  793  443  922  831  866 
1-000000812 

1  ■  000  000  000  793  443  922  831  866 

800()niiG(M)(i.S4  75513S 

10  COO  OOU  007  934  439 

2000000001586  888 

1  -  000  000  812  793  444  567  108  331 
1-000105 

1  •  000  000  812  793  444  567  108  331 

100000  081  279  344  456  71 1 

5  000004  063  967  222  836 

1  -  000 105  812  878  787  878  787  878 
1-320 

1  -  000  105  812  878  787  878  787  878 

300  031  743  863  636  363  636  363 

■20  002  116  257  575  757  575  757 

1  -  3-20  139  672  999  999  999  999  998 


Since  the  characteristic  in  the  given  logarithm  is  1,  there  will  be  2  figures 
before  the  point  in  the  corresponding  antilogarithm.  Thus  the  required 
result  is 

13-  -201  396  7-29  999  999  999  999  98. 

Examjile  J.— Given  7r  =  3"  141  592  653  589  793  238  462  643  show  that 
log  TT  =  0-  497  149  872  694  133  854  351  268. 

Examjile  4.— Show  that  the  antilogarithm  of 

0-434 294 481903-25 18-27  651  129  is  2- 718 '281  8-28  459  045-235  360-290. 

20.  Inverse  Interpolation. 

Suppose  the  values  of  a  function  liavc  been  tabulated,  and  that 
it  is  required  to  find  the  value  of  the  argument  which  corresponds 
to  some  given  value  of  the  function,  intermediate  between  two  of 
the  tabulated  values.  The  process  by  which  we  obtain  this  value, 
which  may  be  called  the  "  autif unction,"  is  known  as  Inverse 
Interpolation. 
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Lct/(a  +  n  w)  be  the  given  value  of  the  function.     Then 

n  (ri  -  1)  -„  , 
/(a  +  n  iv)  =y;  +  n  Sy^  +  -^3  !  '^'  "^  ' '  ' 

r\  '  T       " 

In  this  equation  f{a  +  nw)  is  given,  and  as  the  values  of 
/j,/i,  ...  have  been  tabulated,  the  values  of  5/,,  S-/j,  ...  are  also 
known.  Hence  the  only  unknown  quantity  in  the  equation  is  ?i, 
and  this  is  the  number  required  for  the  determination  of  the  anti- 
function  (a  +  71  w). 

To  obtain  it  we  must  solve  this  equation.  It  has  been  pointed 
out  already  that  the  differences  of  the  function  under  consideration 
usually  become  sensibly  constant  at  some  definite  order,  say  r. 
The  higher  differences  will  therefore  be  zero,  and  this  equation  then 
becomes  one  of  degree  r  in  n.  It  can,  therefore,  be  solved  by  a 
method  of  successive  approximation,  as  follows  : — 
If  we  write  the  above  equation  in  the  form 


«A-^«V> 


and  neglect  differences  of  the  second  and  higher  orders,  we  get  as 
a  first  approximation  the  value  n^  where 

/-/ 

Let  us  next  take  second  differences  into  account  and  put  n^  for  n 
in  the  denominator  of  the  expression  on  the  right-hand  side.  Then 
the  next  approximation  n.,  is  given  by 

n,= l^ . 

When  third  differences  are  taken  into  account,  n.^  is  substituted  in 
the  denominator  for  n-^,  giving 

^  = IzA 

'      S/;-hJ(n,-l)8V;-t-i(n,-l)K-2)SV:' 

Proceeding  in  this  way  we  obtain  closer  and  closer  approxima- 
tions to  the  value  of  n.  The  method,  though  rather  laborious,  has 
the  advantage  that  any  error  introduced  in  the  computation  of  any 
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one  approximation  will  not  cause  the  final  result  to  be  wrong, 
being  itself  rectified  in  the  higher  approximations.  Thus  the  sole 
effect  of  the  error  introduced  is  merely  to  make  the  approximation 
less  rapid,  and  therefore  to  increase  the  amount  of  labour  involved 
in  the  accurate  determination  of  n.  On  this  account  the  process  is 
a  safe  one  to  employ. 

Other  formulae  for  determining  the  antifunction  may  be 
obtained  by  reverting  any  of  the  formulae  of  interpolation.  We 
shall  exemplify  this  by  reverting  Stirling's  formula 

J  (a  +  71  tv)  =/,:,  +  n  /x  S/o  +  ^  d-.t„  +       '  g  ,        M  ^A 

+  r<!r^aY.+  ""''".5V"'''V«v;+... 

Neglecting  differences  of  higher  order  than  the  fifth,  and  arranging 
in  powers  of  n,  we  have 

/(a  +  nw)  =/o  +  (/x  6/o  -  i  /x  SV;  +  ;,V  M  Kfo)  n  +  {\  S'fo  "  -V  ^^0)  n- 

which  may  be  written 

F=An  +  Bri"  +  Cn^  +  Dn''  +  Eiv" 
where 

Reverting  this  series  we  obtain  successively 
F 


F        n_fF\' 
1  ~  ~A 


(i) 

(ii) 

F        B    (F         B  /F\'\'       C   /FY 

<^")  ''=T-  -J  [A-  iXi)}  -:k[a) 


F         B_  (F 
1  ~    1 
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-T'  xU/    rU/    aj [A) 


.  BC 


4    n^/j  ui 


A    more    convenient    form  *    may    be     obtained    by    putting 

F  n 

— -  =  w,  and  -^  =  r  and  rearransfinsr  the  terms. 

A  A 

We  then  obtain  the  formula 

n  =  n,  +/;  nj  +/,  rif  4-/3  ?z/  +/, n/ 
where 

f^=  -Br  +  1  {Brf  -  5  (5r)-^  +  14  {Br)' 
f.  =  Cr{  -  \+b  Br- -21  {Brf} 
/s  =  Dr{-l+6Br)  +  3{Crf 
f\=  -  Er. 

If  the  first  derivate  of  the  function  /(a  +  mv)  is  tabulated,  the 
values  of  the  quantities  A,  B,  C,  D,  E  may  be  obtained  by  forming 
the  successive  differences  of  wf  {a  +  n  w).  As  fourth  diJBferences  of 
the  derivate  are  of  the  same  order  as  fifth  differences  of  the 
function  we  need  not  extend  the  table  beyond  this  order  of 
differences. 

*  Van  Orstrand,  Phil  Mag.  (6)  15  (1908),  p.  630. 
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Argument. 

Entry. 

A 

A-' 

A» 

A^ 

a-3?r 

U\/"_:> 

w5/'_,-. 

a -2  IV 

wr_. 

w  B/'  _ .. 

?r8\f'_., 

.■^f:^ 

a-tv 

"•y'-i 

^oBf_^ 

«-sy'_, 

w5V_'^ 

w8'/'_ 

a 

''vf'n 

10  8/'^ 

w8V', 

v,sy\^ 

TV  8'/'" 

a  +  w 

wf\ 

iv8/\^ 

wo-f\ 

wBV'^ 

w8'f\ 

a  +  2  ?i; 

wf'-i 

7V  8/', 

wB\r, 

a  +  2  w         wf'-i 
Diffei-entiatino;  the  equation 

f{a  +  nw)  =yi  ■{■  An  +  Bn-  +  C'n?  +  Dn^  +  Eiv' 
we  see  that 

A=wf\,;     B  =  \w\f;' ;     C=^lw"f,;";     D  =  ^}fW%^'> 

Now  from  Stirling's  formula  we  have 

wf  {a  +  nw)  =  v)f\  +  nwii8f\+  ~iv8\f\+  — — —  iv  jxBi^f 

n-  (n^  -  1)      -,  ., 
■        ^  'w8'J', 


4! 

o  /-//  /              X            <>  /.,             CO..,       3 «-  -  1         ,,„  ,,       2  «■■  -  n     „.  ., 
7V-f    (a  +  w  ?^•)  =  n-[M8/  „  +  n ?r5;/  ^  +  ^ ?'•  j«'\/  o  H ^o  "  ^'    •'  » 

?6'y"(a  +  n?r)=  ?/;  5-/',,  +       vi  ?/; /x  S''/ 'o      +        n  o"  ""^^'0 

?r*_/'''"  (a  +  n  tt')  =  7^>  /x  8-'\f'„      +  mv  8*/'^ 

?tV""*  (a  +  w  w)  =  iv  8*f\ 

Putting  n  =  0  in  these  equations  we  obtain 

A  =  ivf'{a) 

B  =  \  w-f"  (a)  =  hwix  8/'q  -  ,V  n-  IX  Sy „ 

c  =  1  u^f"  (a)  =  I  tv  sy'o  -  J.  zv  sy '„ 
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A  third  method  of  determining  the  antifunction  is  to  determine 
Wj  as  in  the  first  method,  and  then,  using  shorter  intervals,  to 
retabulate  the  function  for  values  of  n  in  the  neighbourhood  of  n^. 
The  second  approximation  is  now  obtained  from  this  new  table  in 
the  same  way  as  n^  was  obtained  from  the  original  table.  One  or 
two  repetitions  of  this  process  will  give  the  value  of  n  to  a  very 
high  degree  of  accuracy. 

An  approximation  to  the  value  of  the  error  committed  by  computing  n 
by  the  first  method  will  next  be  obtained. 

Since 

/-/o 

n= 

/-/o 


n-\ 


l^i 


the  numerical  value  of  the  error  committed  is 

/-/o      n-l^    „    f"{x  +  n'w)    I 
,     54    •  -2  !     "'"  54  I 

I    n(n-l)        f"  (x  +  n'  w)   I 

Example  1. — The  following  table  gives  the  logarithms  of  the  distance 
of  Venus  from  the  earth.  To  find  when  the  logarithm  had  the  value 
9-9351799. 

A  A2  A3 


Date. 
1914    Aug.  18 

20 


9-9617724 
9-9543645 
9-9468069 
9-9.390950 
9-9312245 
9-9231915 


-  74079 
-75576 

-77119 

(-77912) 

-  78705 

-  80330 

-  81990 


A2 

-1497 
-1543 
-1586 
-1625 
-1660 


-46 

-43 

(-41) 
-39 

-35 


A^ 


A5 


30  9-9149925 

We  shall  use  the  second  method  described  above. 
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Let 

/o  = 

9-9390950. 

Then 

A  = 

-0-0077912  + 

0-0000007=  -0-(^»77!>(»5 

B  = 

-  0  0000793 

C=- 

-0-0000007 

D  = 

0 

E=^ 

0 

11^- 

9-9351799-9-9390950            

-0-0077905          -Ui)ir„o.io 

r  = 

-64-5062 

Br  = 

0  0051153 

Cr  = 

0-000O452 

Dr  = 

0 

Er^ 

0 

fi  = 

-   0-0(D50637 

/,,),  =  -0-00-2545 

fo  = 

-   0-00(X»44O 

./;hi-=  -0-000011 

'A= 

0 

«!=     0-502535 

f'<- 

0 

n  =     0-499979 

Hence  the  time  required  is 

1914     Aug.  (24  d.  +  0-499979  x  43h.) 
or     1914    Aug.  24  d.  23  h.  59  niin.  56  s. 
Example  5.— To  show  that  when  .t=1(X»0  and  w  -c  1  the  error  in  com- 
puting w  from  the  formula 

log  (x  +  w)  -  log  g; 
""""loglx+M-logx 
is  less  than  0-000125. 
Here 

I  ir(l-y)   f"{x  +  w')  I      ^^        ^  ,       , 

I  £r  I  =  , —    — f"^ —     where  0  <  "■  <  I 

I  „-(\-v)         m 1 I 

"^  I         2      '    {x  +  ('-')-    log  (a;  +  1 )  -  log  x  | 

But  by  Taylor's  theorem 

ra 
log(.T+l)-loga-=    ^  _^  ^^,„     where  0-=*"' <1. 

iv(\-)v)      .r  +  w"    I 


\E\ 


{x  +  w'f 


^-     s    •     10  0  0 

<  0-000125. 
Example  5.— The  following  table  gives  the  values  of  y,(.v)  :— 

X  J,(x) 

1-1  0-47090 

1-2  0-498-29 

1-3  0  ■5-2202 

1-4  0-54195 

1  -5  0-55794 

1-6  0-56990 

1-7  0-57777 

Find  the  value  of  x  corresponding  to  the  value  055302  of  Ji{x). 
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21.    Various  Applications  of  Inverse  Interpolation. 

The  method  of  inverse  inteipolation  may  also  be  employed  to 
obtain  the  maximum  and  minimum  values  of  a  function  which  is 
given  as  a  series  of  observations,  or  as  a  graph.  In  this  case  we 
should  differentiate  one  of  the  formulae  of  interpolation,  say, 

f{a  +  n  w)  =M  +  n  S/,  +  ^^^^  Sy,  +  . . . 


Differentiation  with  respect  to  ?^  gives 

2  71  -  I 

wf  (a +  71  w)  =  5 /;  +  — —  Sy  +  . . . 
For  maxima  or  minima  we  have 

/'  {a  +  n  w)  =  0, 
so  that  71  is  given  by  the  equation 

0  =  6/,  +  — -—  S-_/i  + 0-/,  +  ... 

which  is  solved  in  the  same  way  as  the  corresponding  equation  in 
last  section. 

When  it  is  desired  to  obtain  the  value  of  ?i,  for  which  the 
gradient  of  the  function  has  a  particular  value  c,  say,  the  equation 
for  determining  n  becomes 

The  same  method  may  be  employed  to  obtain  the  solution  of 
any  equation  involving  one  unknown,  whether  it  be  algebraic  or 
transcendental.  For  we  may  tabulate  the  values  of  the  expression 
occurring  in  the  equation  and  obtain  their  successive  differences. 
The  solution  of  the  equation  is  then  the  same  as  finding  the  values 
of  the  argument  corresponding  to  zero  values  of  the  function. 

Example  i.— Approximate  to  the  roots  of  the  equation 
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First  form  the  table. 


X 

/(-r) 

A 

A2 

A' 

-3 

-   2 

1-875 

-2-5 

-   0125 

01-25 

-  1  -750 

0-75 

-2 

0 

-   0-875 

-  1-000 

0-75 

-1-5 

-   0-875 

-    1125 

-0-250 

0-75 

-1 

_   2 

-   0-625 

0-500 

0-75 

-0-5 

-    2-625 

0-625 

1  -250 

0-75 

0 

_    o 

2-6-25 

2-ftOO 

0-75 

0-5 

O-G-25 

5-375 

2-750 

0-75 

1-0 

6 

8-875  ♦ 

3-500 

0-75 

1-5 

14-875 

13-125 

4-250 

0-75 

2  0 

28 

18-125 

5-000 

0-75 

2-5 

46  125 

23-875 

5-750 

3  0 

70 

Tlie  positions  of  the  roots  will  be  found  by  examining  the  column  J'{.r). 
We  then  see  that  -when  /(:r)=:0,  the  value  of  .c  lies  between  0  and  05.  It 
will  also  be  seen  that  f{x)  is  zero  when  x-  -2,  and  that  there  is  a 
possibility  of  another  root  in  the  interval  -2-5<  .r  <  -  1-5.  We  shall 
consider  the  various  approximations  to  the  root  for  which  0<  a-  <  0  5. 

Using  the  formulae  of  the  first  method  of  §  20  we  have 

0-(-2) 


- 

2-6-25 

= 

0-7619 

/-/o 

~ 

5/j  +  i(ni- 

DS-Vi 

2 

-2-625- 0119  X  2-75 

=  0-8704 


60  INTERPOLATION  AND  NUMERICAL  INTEGRATION    [cii.  in 

f^ 

^■^  ~  sj\ + M"2  - 1 )  5v\ + H  »i2  - 1 )  ( «2  -  2)  5v; 

2 

^  2-625  -  0  0648  x  2  75  +  0  0244  x  0  75 

=  0-8113 


'  "  5A  +  i(H;i-  1)  dy]+h{",-l)  (n,-2)  5V, 

"  2^625  -00944  x  2-75  +  0-0374  x  0-75 
=  0-8356. 

The  first  four  approximations  to  the  value  of  n  are  therefore 
?«!  =  0-7619 
no  =  0-8704 

713  =  0-8113 

714  =  0-8356 

But  the  tabular  interval  is  one-half.     Hence  the  first  four  approximations  to 
this  root  are 

Xi  =  0-3810 

X2  =  0-4352 

xg  =  0-4057 

iv,  =  0-4178 

The  correct   value  of   the   root   is  0-414-2.     Subtracting  each  of  the  values 
obtained  from  04142  we  get 

0-4142 -.T]  =  +0-0332 
0-4142 -a'.2  =  -0-0210 
0-4142 -3-3=  +0-0085 
0-4142-.T4  =   -0-0036 

These  differences  show  that  each  step  gives  a  much  closer  approximation 
to  the  root,  and  that  by  continuing  the  process  we  can  get  a  result  which 
shall  be  correct  to  any  given  number  of  decimal  places. 

Examples  of  the  second  method,  in  which  only  first  differences  are  taken 
into  account,  will  be  found  in  any  text-book  which  discusses  the  solution  of 
numerical  equations. 

Example  2. — Approximate  to  the  real  roots  of  the  equations 
(i)  e-^-.T3=:0 

(ii)  e^  +  a;2-4  =  0 

(iii)        0-5  x-i-5  -  12  logio  a'  +  2  sin  2  a:  =  0-921. 
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MISCELLANEOUS  EXAMPLES. 

1.  If  ^  =  .1- tan-^  .r,  calculate  a  table  of  first  differences  near  .i-  =  0-4.3  for 
increments  00002  in  x: 

2.  The  following  table  gives  the  moon's  right  ascension 


Date. 

R 

A. 

,  March  5d. 

(Ih. 

5h.     :]m. 

7  "26  8. 

Kill. 

26  m. 

58-96  s. 

20  h. 

51  m. 

14-83S. 

6d. 

6h. 

6h.    1.5  m. 

48-96  8. 

10  h. 

4(1  m. 

.34 -76  s. 

7(1. 

2h. 

7h.     5  m. 

25-36  s. 

Complete  the  table  for  every  hour  from  March  5(1.  20  h.  to  March  6d.  6  h. 

.3.  Prove  that  in  a  table  of  logarithmic  tangents  to  base  10,  the  difference 
for  one  minute  in  the  neighbourhood  of  60^  is  0-00029. 

4.  Determine  the  values  of  (1)  log  tan  1    29'  3.S" 

(2)  log  sin  V  15'  12" 

5.  Show  that  logi,,  1031  =  3013  258  665  -283  516  546  909  664. 

6.  Show  that  the  number  whose  logarithm  to  20  places  of  decimals  is 
given  by  0-  004  892  890  303  -239  039  78  is  1  -01 133. 

7.  From  the  following  table  determine  the  value  of  x  fur  which  /(.'.)  has 
the  value  0-22.389. 

;t-  /U) 

1-5  0-51183 

1-7  0-39798 

1-9  0-28182 

2-1  0-16661 

2-3  0-05554 

8.  Appro.ximatc  to  the  real  roots  of  the  equations 


(i) 

x-e^  +  2.i--5  =  0. 

(ii) 

e'-e— '  +  0-4.t;-  10-0 

(iii) 

2.i-3-i-3.t--16-0 

(iv) 

2-42.r- 3-15  log^..c -20-5-0. 

From  the  series 

berx-  =  1 


2-' .  4'-^  ^  2^ .  4- .  6- .  8-  ~  2- .  4'-^ .  6- .  8'^ .  10- .  12- 


compute  a  table  of  ber  a:  from  x-  =  0-0  to  a,-  =  5-0  at  intervals  of  05,  correctly 
to  6  places  of  decimals. 

Hence,  by  subtabulation,  compute  a  table  of  ber  .f  from  a.-  =  2-0  to  x  =  3-0 
at  intervals  of  O'l. 
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10.  The  confluent  hypergeometric  function  W^,  ^^^(z)  is  computed  for  large 
values  of  z  from  its  asymptotic  expansion 

while  for  small  values  of  z  it  is  computed  from  the  formula 

r(-2m)  r(2m) 

^., .. (^)  =  Til-m-k)  ^^  .. ^^)  +  r(Km-)^)  ^^.  -... (^) 

where 

(  h  +  m-  k  (^  +  m  -  k)  {j  +  m-  k)  \ 

M,_  ,„  (z)  =  z^+'^  e-5^|l  +  li  (2m  +  l)  ~  +  2!  (2m  +  lH2m  +  2)  ^■"+  -j  _ 

Compute  tables  and  draw  graphs  of   W^.  ^Jz)  for  positive  values  of  z  in  the 

^•=-3-1,  7n=+2-2 
k=-0-l,  7n=+()-2 
k^+2-2,     711= +0-2. 


CHAPTER    IV 

NUMERICAL   INTEGRATION 

22.  Introduction. 

The  method  of  evahiatiiig  the  delinite  integral  of  a  function 
from  a  series  of  numerical  values  of  the  function  is  called 
Mechanical  Quadrature  or  Numerical  Integration. 

If  a  function  is  given  by  its  graph,  the  simplest  way  of  deter- 
mining an  area  bounded  by  the  curve,  two  given  ordinates  and  a 
given  abscissa,  is  by  means  of  a  planimeter.  This  method  would 
naturally  be  used  in  the  case  of  the  indicator  diagrams  of  steam, 
gas,  or  oil  engines,  or  the  cards  from  certain  hydraulic  motors, 
or  the  stress-strain  diagrams  drawn  by  various  types  of  testing 
machines.  The  accuracy  of  the  result  obtained  would  in  such  ^ 
cases  be  as  good  as  is  required,  but  it  is  to  be  remembered  that 
graphical  methods  are  not  susceptible  of  great  refinement.  If, 
therefore,  considerable  accuracy  is  required,  or  indeed  in  general 
whenever  the  function  is  specified  by  a  table  of  numerical  values,  the 
method  of  numerical  integration  is  preferable  to  the  use  of  the  plani- 
meter. The  method  of  numerical  integration  must  also  be  resorted  to 
when  a  function  is  specified  by  its  analytical  expression,  but  cannot 
be  integrated  in  terms  of  known  functions  by  the  methods  of 
the  Integral  Calculus  ;  e.g.,  the  elliptic  integrals  belong  to  the  class 
of  functions  which  cannot  be  evaluated  by  the  elementary 
methods  of  the  integral  calculus. 

The  formulae  necessary  for  numerical  integration  are  derived 
from  the  formulae  already  established  for  interpolation.  As  in 
the  case  of  interpolation,  the  order  of  differences  which  must  be 
taken  into  account,  will  depend  entirely  on  the  rapidity  with 
which  the  differences  decrease  as  the  order  increases.  It  need 
hardly  be  said  that  unless  the  convergence  of  the  series  can  be 
ensured   the  process  is  of  no  avail.     Consider,   for  example,   the 
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determination  of  the  area  of  a  semicircle,  the  ordinates  being 
perpendicular  to  the  diameter.  No  matter  how  many  ordinates 
are  taken  (i.e.  no  matter  the  order  of  differences  employed),  the 
approximation  to  the  area  by  means  of  many  of  the  formulae  to  be 
hereafter  established  will  not  be  a  good  one.  But  if  we  evaluate 
the  area  included  between  an  arc  of  a  circle,  less  than  a  semicircle, 
a  diameter,  and  the  two  perpendiculars  from  the  extremities  of  the 
arc  on  this  diameter,  we  can,  from  the  same  formulae,  obtain  an 
approximation  to  the  true  value  with  any  degree  of  accuracy  we 
please.  The  reason  for  the  failure  in  the  first  instance  is  that  the 
differential  coefficients  of  the  ordinates  at  the  limits  of  integration 
are  infinite,  and  as  a  result  the  series  is  not  convergent.  For  this 
same  reason  we  must  assume  that  the  function  represented  by  the 
observations,  and  its  differential  coefficients  up  to  the  order  of 
differences  employed,  are  continuous ;  a  condition  which  is  not 
necessarily  fulfilled  in  natural  problems,  as  may  be  seen  by  an 
examination  of  an  ordinary  barometric  curve. 

23.  Evaluation  of  Integrals  by  the  Integration  of 
Infinite  Series  Term  by  Term. 

Before  applying  the  methods  of  finite  differences  we  shall  give 
an  example  showing  how  in  many  cases  the  value  of  the  integral 
of  a  function  which  is  analytically  known  may  be  obtained  by 
expanding  the  function  in  a  uniformly  convergent  infinite  series, 
and  then  integrating  it  term  by  term. 

Example  1. — Consider  the  complete  elliptic  integral 


^^-^<W^i 


sin"  X ) 


=  I    dx(\  -k^  sin-  x)   ^ 

Jo 


Assuming  that    |  A;  |  <  1  and  expanding  the  term  (1  -  k-  sin"  x) 
by  the  binomial  theorem,  we  have 

IT 

^  (^>  ^)  =  [ "c^  a;|l  +  I  yi-  sin^  x  +  ,^— ^  k'  sin^  x 

Jo         '^         "  Ij".  w  . 

+  k^  sm*^  X  +...|  . 


-i 
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1  .  3.5  ...  (2n-l)  TT 


But 
Hence 


If  k  is  .sinali  tliis  .series  is  very  suitable  for  obtaining  the 
\alu(>  of  F  ;  but  it  is  evident  that  if  k  appi'oaches  unity  a  hiri,'e 
uumboi-  of  terms  must  be  included  to  i;et  a  fair  approxima- 
tion to  tlie  value  of  /•'. 

Let  ^■  =  0'1.  Then  to  obtain  F  correct  to  4  decimal  places  it  i.«  only  neces- 
sary to  include  two  terms  of  the  series.     For 


/  ,r\       TT  )        0--J       0140625 


100  ^    10000    ^■■- ) 

=  ^{1+00025 +  0  0000 140625  +  ...} 
=  ^{1-0025} 
=  1-5747. 

Example  i".— Find  the  value  of  F(^I:,  ^)  when  i-  =  sin  10° 

IT 

Examplf:  3. — Show  that  I  ^X-IS-^m-x  is  equal  to  r")44l5when  k  has 
the  value  sin  15°. 

2     ro-6 

Example   4. — Show   that    the  value   of  the   integral  — 7^^   I       t~'^ dx  is 

V  'T    Jo 
0  60386. 

24.  A  Formula  of  Integration  Based  on  Newton's 
Formula. 

We  shall  now  show  how  integration  may  be  performed  by  tlic 
aid  of  interpolation-formulae.  Intcgiating  Newton's  interjiolatioii 
formula 

J{a  +  n  w)  =/o  +  n  c>/^  +  ^^-^-^ — -^  ^-J\  +  . . . 
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with  respect  to  n,  we  have 

\f{a  +  7i.n:)dn  =  Jl  \dn  +  8/,  \ndn+  ^  ^Vi  \n{n-l)dn 

+  ^^8yAu{7i-l){7i-2)dn+... 

Hence 

f /(«  +  n  w)  dn  = /;,  +  h  S /;  -^^J\f,  +  h  ^A  -  T^^V  s'./2  +  •  •  • 
Jo 

f /(«  +  n  w)  dn  = /;  +  h  5. 4  -  tV  ^Vi  +  zh  ^'/#  -  t^  ^'/^  +  •  •  • 

f /(«  +  «  w)  dn = /:  + 1 8.4  -  Jj  sv;  +  ^  svj  -  T^^V  s^/  +  •  •  • 
["/(«  +  7i «')  d  n = /;„,  + 1 8 /;_:  -  T^,  sv; + ^  sv;+^  -  ,w  s^y:.+i  + . . . 

Jr-l 

Adding  these  equations  we  have 
f{a  +  ri  ?f)  c?7i  =    /o  +/,  +  . . .  +./;_! 

+  MS4+S/3+  ...  +sy:_,) 

-TV(syi+sV.+ ...  +sv;) 
+  .i^(sVn+8=i/;+ ...  +s^./:.^i) 


But 

syi  +  8/.+  ...  +  8/,_,  =/;  -y;  +./;  -./;  +  •  •  •  +./;■  -/-I 
=/:.-/o 

sv» + s^/. + . . . + 8v;.+i = sy,  -  sy, + s%  -  sv;  +  •  •  • + sy+1  -  sy 

'  =  sy^i-sy 
sy+sv3+...  +  sy+i  =  s!/^-8!/;+s^A-8!/k+ ...  +%.,-s3/;_^. 


^r 
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Hence 

/{a  + 11  ?(■)  dn  =  Ay; + /;  + /;  + . . .  +/,_,  +  iy; 

-i'M^'/r+^-^/^+ 

Putting  x^a  +  nu-  this  gives 

/{x)  dx  =  i/i  +./;  + /;  + . . .  + /;_,  +  y; 

-rVU57;.+«-5Vi)  + 

This  expansion  will  give  the  value  of  the  definite  integral,  pro- 
vided it  is  possible  to  obtain  values  of  f(x)  outside  the  limits  of 
integration  in  order  to  obtain  the  differences. 

When  this  result  is  written  in  the  foim 

^1  ■/(•«)o?«^=(./;+/;+---+y;_,)+i(./;+/o)-A(^./;4-i-5/;) 
+ ^  (S7;+i  -  ^vi)  -  ^  (sv;.+e  -  ^v;) + •  •  • 

it  is  easily  seen  that  the  coetJicients  are  those  of  x  in  the  expansion 
1  +  A(,  X  +  ^1  X-  +  ^,,  a;^  +  . . . 


I 


log(l+.T) 

Clausen  *  has  computed  the  first  thirteen  of  these,  v; 

'~  "^T  ''  ~  ~  3628800 

1  ,  8183 


A. 


12  °    1036800 

I  .     3250433 


24  "  479001600 


19  4671 

^z~  ~  ;rT7,  -"10 


720  '"    788480 

_3_  13695779093 

'""^160  "~  ~2fil534S736000 

863  ,      2224234463 

Ar.=    -  -r-——:  il.o 


^. 


60480  '-    475517952000 

275 


2419: 


Jour,  reiiie  angeic.  Math.,  6  (1830),  p.  287-9, 
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In  the  same  memoir  it  is  shown  that  when  the  limits  of 
integration  are  a-^w  and  a  +  {r  +  h)  w  the  coefficients  are  those 
in  the  expansion 

— ^ =  1  +  B^x-  +  B..af'  +  B.,x^  +  ... 

J{\+x)\og{\+x) 

These  are  also  tabulated. 

The  formula  obtained  above  must  be  subjected  to  another 
transformation  when  values  outside  the  limits  of  integration  cannot 
be  obtained. 

Now 

(i)  s./;+,  =  sv:.  +  sy;.-. 

(ii)  sy,+i-2S7;  +  sv;._,  =  sv:. 

i.e.        sy,+,-2sy,=  -s\t]._,  +  ^^f, 

(iii)  st;.  =  s^y;^,  -  3  sv;+,  +  3  s^y;,,  -  8v;_ , 

Substituting  these  values  in  the  formula  we  have 


tVja 


fix)  dx  =  \f,  +/  +/,  + . . .  +y._,  +  ^/. 


This  form  involves  only  differences  which  can  be  calculated 
from  values  of  the  function  which  are  within  the  limits  of 
integration,  and  is  the  form  to  be  used  when  no  values  outside 
those  limits  can  be  obtained. 

Example  i.— To  determine  w  from  the  equation 
TT         r  ^      dx 

using  differences  of  J^. 

Here  f{x)  =  j-^j-^  ;    f,  =/(0)  =  1  ;    /,=/(!)  =  0-5  ; 

i«  =  01  ;     r  =  10. 

Inspection  of  the  formula  will  show  that  it  is  not  necessary  to  form  a 
complete  table  of  differences.  The  highest  order  of  differences  retained  will 
be  the  4th,  and  as  5^/2  ^'^'^  ^^^8  ^"^s  the  only  two  of  this  order  required,  the 
central  portion  of  the  table  may  be  omitted  as  in  the  following  scheme  : — 
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X 

/W 

A 

0-0 

1-00000 

-  990 

0-1 

0-99010 

-2856 

0-2 

0-96154 

-4411 

0-3 

0-91743 

-5536 

0-4 

0-86207 

0-5 

0-80000 

0-6 

U -73529 

-6415 

0-7 

0-67114 

-6138 

0-8 

0-60976 

-5727 

0-9 

0  •55-249 

-5249 

1-0 

0-50000 

From  this  table  we  have 

i/u  =  0-50000 

J\  ^  0-99010 

/_,  =  0-96154 

/,  =  0-91743 

/,  ==  0-86207 

/,  =  0-80000 

/,  =  0-73529 

/-  =  0-67114 

/,  =  0-60976 

/g  =  0-55249 

i/„=  0  25000 

.Sum  =  7-84982 

Hence 


-1866 
—  1555 
-11-25 


A=» 


311 
430 


•277 
411 

478 


134 
67 


-0  05249 
-0-00990 

Diff.  = 

-0-04259 

0  00478 
-0-01866 

Sum  = 

-0-01.388 

Diff.  = 

0  00067 
0  00311 

-0  0(.)244 

5Vh    = 
5'/.    = 

-0-00067 
0  001 19 

Sum  = 

0  00052 

'1      dx 


•84982   -    -^  (-0-04259)    - 


(-0-01388) 


A* 


119 


160 


(0-t  10052) 
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=  7-84982 

355 

58 

6   -    I 

=  7-85400 

IOtt 

—r-  =  7-85400 
4 

TT  =  3-14160. 

rs      xdx 
Example  ;?.— Determine  loge  2   from  the  equation  h  log^  2  ^  l        ^  ^  ^.2 

using  differences  of  ^'jj. 

25.  Case  when  the  Upper  Limit  is  not  a  Tabulated 
Value. 

In  the  previous  section  it  has  been  assumed  that  the  quantity  r 
in  the  upper  limit  was  an  integer,  so  that  the  tabular  value  of  f{x) 
was  known  when  x  had  the  value  a  +  riv.  When  the  upper  limit 
is  a  +  r'w,  where  r<?-'<r+l,  and  the  lower  limit  coincides  with 
one  of  the  values  of  the  argument,  so  that /(a)  is  known  when 
x  =  a,  a  +  'u;  ...  a  +  riv,  a  +  r  +  lic,  ...  the  value  of  the  integral 
may  be  obtained  as  follows  : — 

First  determine,  as  in  last  section,  the  value  of 

f{x)dx. 

Then,  dividing  the  range  {a-\-rn\  a  +  r'ir)  into  a  convenient 
number  of  equal  intervals,  interpolate  the  values  of  J'{x)  for  each 
of  these  sub-intervals.     The  value  of 

f  "+'■«' 

f{x)  dx 

Ja  +  rw 

may  now  be  obtained  by  making  use  of  these  new  values  of  f(x) 
and  their  differences. 

When  the  lower  limit  does  not  coincide  with  a  tabular  value  of 
the  argument,  a  similar  process  will  enable  the  value  of  the  integral 
between  the  given  limit  and  the  nearest  tabular  value  of  the 
argument  to  be  obtained. 
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Another  method  which  might  be  employed  is  to  obtain  the 
value  of  I  f{x)dx  for  different  ranges  of  integration,  and  then 
determine  the  value  for  the  given  range  by  interpolation. 

CZV  20' 
Example  1. — Evaluate  I  aosxdx,   given   the  values  of  cos  x  when 

J -20° 
x  =  20°,  22°,  24",  -26%  28%  30%  32°. 

Forming  the  table  of  differences  we  have 


X 

cosx- 

A 

A- 

A' 

20° 

0-9396926 

-  125087 

2-2° 

0-9271839 

- 136384 

-11-297 

166 

24° 

0-9135455 

- 147515 

-11131 

182 

26° 

0-8987940 

- 158464 

-  10949 

191 

28° 

0-88-29476 

-  169222 

-  10758 

2U7 

30° 

0-8660254 

- 179773 

-  10551 

32 

0-8480481 

Fir>it  Method. — We  shall  first  of  all  determine  I       cos  xdx.     Substituting 

J  20" 
in  the  formula 

i-  f '   '  '\f{x)  d  X  -  i/o  +/i  +  •  •  ■  +/.-1  +  \fr  -  h  (S/, -  .  -  5y\ )  -  .h  (5-/,-i  +  5Vi) 


we  have 
0 


1  P^"  /       n.^fioo^Ao      .  r -169222-',      ,  r-lOTDSn       i ,,   f    mi 

034M66j,o.°°^^-'^'^  =  nff,f-'-'L 


0- 

0-9271839 
0-9135455 
0-8987940 
0-8829476 
0-4330127 

=  4-5-257896 

50" 
cosa;cZa; 

20" 

=  01579798. 

Now  by  using  Newton's  formula  for  backward  interpolation 

n{n-\)                n(n-\)(n-2)    „ 
f(a-n  „')  =/„  -  n  d/_  ,  +      \^  y  '  5 V- 1  "  ^ gT S'V  -  3  +  - 

determine  the  values  of /(.i),  i.e.  cos  x,  when  .t-  has  the  values  30°  -20',  30°  40' 
3r,  3r  -20',  31°  40'.     We  thus  obtain  the  following  table  :— 


INTERPOLATION  AND  NUMERICAL  INTEGRATION     [cu. 


X 

SCO' 
30'  20' 
30'  40' 

sro' 

3r20' 

;ir  40' 

32°  0' 


cos.r 
0-8660254 

0-8631019 
0-8601491 
0-S571673 
0-8541564 
0-85II167 
0-8480481 


-  29235 

-  295-28 
-29818 
-30109 

-  30397 

-  30686 


-  -293 

-290 
-291 

-288 
-289 


From  this  table  we  obtain 
1 


0  0058178 


f  ''"  ^os  .c  cZ.^0-43301-27  -  .J  "  ^^'^'^l  -  -.■  \  "  ''^'^ 
Jao  o.,g3,o^,      -L^-29235j      ^[-293] 


0-8601491 

0-8571673 

0-4270782 

=  3-4405189 

SI''  20' 

cosxrfa;  =  0-0200161 


Adding  these  results  together  we  get 


31"  20' 

cos  ,(■(/.(■  =  0-1' 


959 


Second  Method.  — Oht&m  by  the  general  method  the  values  of   I    cos  .r  (Z.i- 

mit 

i  th 

r 


when  the  lower  limit  is  20"  and  the  upper  limits  are  2-2%  24',  26",  28°,  30%  3-2'' 
respectively.     We  thus  obtain  the  following  table  :— 


20° 
22° 
24" 
26° 
28° 
30° 
32° 


cos  X  dx 
0  0000000 
0-0325865 
0-0647165 
0-0963510 
0-1274515 
0-1579799 
0-1878992 


325865 
321300 
316345 
311005 
305284 
299193 


-  4565 

-  4955 
-5340 
-5721 
-6091 


-390 
-385 
-381 
-370 
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Using  Newton's  formula  of  backward  interpolation  we  obtain  the  value  of 
the  required  integral,  viz. 

cos.cdi-  =  0-187899-J-i(-299193)  +  H7T-  (-6091)-  .^  ,    '(-370) 

=  01779961. 

The  values  obtained  by  the  two  methods  agree  to  6  places  of  decimals. 
The  result  correct  to  S  places  is  0-17799598,  but  as  the  last  figure  is  always 
forced,  the  one  method  can  hardly  be  assumed  better  than  the  other.  Probably 
the  latter  method  is  to  be  preferred,  as  in  the  former  the  division  of  the  range 
between  the  given  upper  limit  and  the  new  upper  limit  into  suitable  intervals 
niuy  give  rise  to  values  of  the  argument  consisting  of  manj*  decimal  figures 
which  would  increase  the  amount  of  labour  involved  in  the  interpolation. 

I*  3U"  -JO' 

Example  2. — Evaluate    I  sina-(/.c,  given    the   values    of   siu.c    when 

J  20« 
x-  =  20°,  2-2°,  24°,  26°,  28',  30°,  32°. 

f-i-5    dx  1 

AVamp/e  J.— Evaluate    I  ,      if  the  values  of         , . .    are   given  for 

integral  values  of  x. 


26.  A  Formula  Based  on  Bessel's  Expansion. 

Anutlier  serie.s  for  the  deliiiite  integral  is  obtained  liy  integrating 
Bessel'.s  formula  of  interpolation. 

"{"-')      >,  .       n{n- \){n  -  }y)  ..,  , 
/{a  +  n  w)  =  ixj,  +  (n  -  h)  S/,  +      '^,        /x  67,  +  -^ -^ &'/^ 

Putting  x  =  a  +  7iiv  we  have  by  integration 

10  ]n  Jo' 

fi  f  1  f  11(11  -  1 ) 

=.  /x /■  \    dn  +  ^./^     {n-  i )  d  n  +  /x  5-./ j  .-,  ,       dn 

'Jo  Jq  Jo       -'  • 


71  (11  -   I)  (it  -   .',)     , 

"-  (/  It 


,^  .fi/i(/r  -  l)(w-2) 


'^i T 


du  + 


/^yi-iV/^«'./i+7yo/^sVi+- 
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Similarly, 

J  «+» 

—  /{^)  dx  =  i^/r-i-  yV  /x  s-y,  _ .  +  j\}^  jx  sv;  _  i  + . . . 

Hence 

—  I  /(x)dx  =  fiJ^  +  fM/,+  ...+fxf^_^ 

+ 

+  tVo  (i  sv;  +  sv;  +  5v:  + . . .  +  av;._,  +  h  svi-} 
+ 

or,  since 

J  2  ^  »  -  a  J^  -  5 

and  therefore 
we  have 

^\7  '"'^'^*^  '^  ^^  ^  --^^  "^-^^  ^-^^ + •  •  •  ^'^'-  ^  ^z-- 

+ 

But 

=  -H^A.  +  sy;) 
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and,  similarly, 

Hence,   finally,  we  obtain  the  formula 

1  f" +'■"■ 
—  /(j:)  d  X  =  hA  +J\  +/■!  +  ■  •  •  +. /i-i  +  hfr 


75 


flO-3  r/.r 

Example  1. — Calculate     I        —  to  8  places  of  d^ 
J  100  •'- 

1 


ecimals. 


Here  f(x) 


/o=100,  r  =  a,  w^\. 


In  order  to  obtain  the  dififerences  required,  it  will  be  necessary  to  form 
f(x)  for  all  integral  values  of  x  from  x-  =  98  to  .c:^107.  These  may  be 
obtained  from  Barlow's  Tables.     We  thus  obtain  the  following  : — 


X 

f(x) 

A 

A- 

A3 

98 

0-010204082 

1(13072 

99 

0(U()101(iln 

-  KilOKl 

2062 

-62 

100 

()'(fl(l(IIKI<HHI 

{  -  KM  1010) 
-  99010 

'2(XK) 

(  -  CO) 
-58 

101 

0  0O99O(J99U 

-  97068 

1942 

-58 

102 

0-009803922 

-  95184 

1884 

-53 

103 

0-009708738 

-  93353 

1831 

-53 

104 

0  009615385 

-  91575 

1778 

-  51 

105 

()-(M»y5-238H) 

(-  90712) 
-  89848 

1727 

(-49) 
47 

lOG 

0(109433962 

-    88168 

1680 

107 

0  009345794 

Substituting  in  the  formula 
'105  iJx 

~-  ^  L'"  ■  A  <  A  ^  A::    j\  i  hA-.  ■  i\  (M  5./;v  -  At  5./;)  +  t",',,  (m  oVis  -  M  s^yi) 


J  100 
we  get 


■  "Ja  dz 

-^0048790940 -0-(XHKXMI775     OOiMyMM^-MMI 
100    ^ 


-0-048790165. 

The  value  correct  to  lo  ijlaces  is  (JO4S79O104-. 
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iVb<e.— It  is  to  be  remarked  that  the  expression 

i/i  +./'i  +/•!  +L  +/.  +  hA,  =  0-048790940 
gives  the  answer  correct  to  six  decimal  places  ;  this  approximation  is  quite 
sufficient  for  most  practical  purposes. 

f '-^  dx 

Example  2. — Approximate  to  the  value  of  the  integral  I        //  .3_q  "ios' 

27.  The  Euler-Maclaurin  Expansion. 

In  the  foregoing  the  value  of  the  integral  is  expressed  as  a 
series  involving  differences.  In  many  cases  in  which  the  function 
is  known  analytically,  it  will  be  found  more  convenient  to  use  what 
is  known  as  the  Euler-Maclauriii  expansion  which  employs 
differential  coefficients  instead  of  differences.  The  rapidity  of 
convergence  is  of  the  same  order  as  that  given  by  Bessel's  formula. 
It  may  be  obtained  as  follows  : — 

Since 

=  h  {J{a+u))-f{a-w)} 

=  l{/{a)  +  ivf'{a)+  |^/"(a.)  +  ... 

-/(a)  +  w/'(a)-  ^/"  («)  +  ...} 

=  w/'{a)+  ^/"'(a)+... 
and  similarly 

/x  S'V;  =  tvy'"  (a)  +  ~  f^"^  («)+... 

jx  5/,  =  wj '  {a  +  rw)+  —f '"  (a  +  r  ?<;)  +  .. . 

fx  f^-'/.  =  I,//'"  {a  +  ?•  w)  +  —  /i^-'  (a  +  riv)  +  ... 
we  have  by  substituting  in  the  formula 

1  f  «+'•'"    . 

— J  ^^     J  (x)  dx^  1 /;  +/  +./;  + . . .  +/._!  +  \j^ 
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the  result 

-  A  ^'•(./;'-/;')+Tk^^"  (/'"-./;'") 

This  is  the  EulerMadaurin*  expansion.  'J'lic  L-oetlicicnts  of  tlie 
various  items  in  it  are  proportional  to  the  well-known  13ernoullian 
numbers.  For  iff  we  put /(x)  =  e'-"  in  the  formula  and  at  the 
same  time  change  the  limits  to  a  -iv  and  a  +  w  ^ve  get 


!  constants. 

—  I         e"-"  dx  =  —  I 


where  A,  7>,  ...  are  constants 
But 


dx 
Hence 

—  («"■  - «-"'.)  =  i/o  +./;  +  hf-i + ^^  "•  (./i'  -.AV)  +  />  ''•'  (./;'"  -/<:")  ^ 

=  ^  e-"  +  1  +  ^  e"'  +  ^  ?'•  (e'"  -  e-'";  +  i?  ir''  {e'"  -  e'") 
from  which  it  is  easy  to  deduce  that 

-z z. -,  =  I  -  Alt"-  ~  B  /(■••  -  . . . 


— -.  cot  7— .  =   \  -  A  7tr  -  B  w*  -  ... 

But  the  BernouUian  numbers  are  the  quantities  B^.  7/.,  ...  in  the 
expansion 

and  therefore  A,  B,  ...  are  proportional  to  B^,  B^,  ... 


*  For    a   rigorous    discussion,    rf.     Whittaker    and    Watson's    Modern 
Analysis,  p.  128. 

t  Poisson,  M6m.  de  I' Acad,  des  Sc,  L^23. 


195 


INTKRPOL.VnON  AND  NUMKRICAL  INTKaRATION     [c...  iv 

Calculate   1        —  to  8  places  of  decimals,  and   check  the 
J  95      ^ 
result  by  computing  it  as  -  log.(l  -  t^)  by  the  logarithmic  series. 

Here  f(,  =  ^,   ^=1,  r  =  5. 
Then 

n^dx_i     1  +i_  +  i.+J_  +  l.+_L.± 

T'gs"  96   97   98   99   2   100 

-  tV  (  -  IW^  +  9p)  ^  ^"^  "  (  "  TTio^  ^  95"V 

^  /  -  0-000100000\  1  _/  -  0000000010\ 
0-005263158-  Tir^ +  0-000110803/  ^  -  "V -i  0  000000012/ 
0-010116667 
0-010309278 
0-010204082 
0-010101010 
0-005000000 
0-051294195 
-  0-000000900 
0-051293295 


Again 


(0-05/   (0-05)3  (0-05)4   (0-05)5   (o-05)6 


-l«g^(l-100r   005  +  -^-+-^ +-^+-^4-   ^ 

=     0-050000000 

0-001250000 

0-000041667 

0-000001563 

0-000000063 

0-000000003 

=     0-051293296 

The  two  results  will  be  seen  to  agree  to  8  places  of  decimals. 

fl05  dx 
Example  2. — Evaluate   I       —  to  8  places  of  decimals. 

J  100    '^^ 

ft 

Example  3. — Approximate  to  the  value  of   i      cos'xdx. 

28.    An  Exceptional  Case- 

Legendre*  remarked  that  if  the  odd  differential  coefficients  above  a 
certain  order  take  the  same  value  at  both  limits,  the  formula  fails  to  give  an 
accurate  value  of  the  definite  integral.     For  instance,  if  the  integral  under 

*  Fonctions  Elliptiques,  Vol.  ii.,  p.  57. 
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discussion  were  the  elliptic  integral  of  the  second  kind  taken  between  limits 
0  and  7r/2,  viz. 


I        J  1-Fsin-  xdx, 
Jo 

all  the  odd  differential  coefficients  become  zero  at  both  limits,  and  thus  the 
value  of  the  integral  would  appear  to  be  given  by 

But  the  value  given  by  this  expression  is  not  in  close  agreement  with  the 
true  value  of  the  integral.  The  reason  of  this  is  that  the  numerical 
coelKcieiits  introduced  in  the  successive  differentiations  increase  without 
limit,  so  that  each  term  takes  the  form  cc  x  0,  which  is,  of  course,  inde- 
terminate. 

Better  methods  of  evaluating  the  three  kinds  of  elliptic  integrals  will  be 
given  in  another  tract  in  the  present  series. 

29.    The  Formulae  of  Woolhouse  and  Lubbock 

We  shall  next  derive  from  the  Eiiler-Maclaiirin  expan.sion  a 
formula  which  is  of  u.se  when  it  is  required  to  evaluate  the  sum  of 
a  large  number  of  terms. 

If  the  interval  between  consecutive  values  of  the  argument  in 
the  Euler-Maclaurin  formula  is  w,  we  have 

—  \y{^)  dx=j\  +/ + /;  + . . .  +./;  - ,',  (./;  +/.)  -  -,v.  v  (/;  -/;) 

W  J  a 

and  if  this  interval  is  ir /  m  the  formula  becomes 

-  [ """' V(^)  dX=/o+/^+f-^+...  + /;  -  h  (/o  +fr) 
'''  J  a  VI  m 

-A^<x'-./i')+.k5'./r -./:")  +  ... 

Subtracting  the  latter  from  ni  times  the  former  we  deduce 

ni  -  I 


ni-  -  1    tr 
12     ^'^' 

-./;■) 

"''-'•'^,v 

■  -  c 

720     m-' 
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This  is  Woolhouse's  Formula.  Obviously  b\"  using  the  right- 
hand  expression  instead  of  the  left-hand  one,  the  labour  involved  in 
determining  the  sum  /^ 4-/2  +./  ^  +  ■••  +/y-  ^^  enormously  reduced. 

A  somewhat  similar  formula,  in  which  differences  are  employed 
instead  of  differential  coefficients,  is  due  to  Lubbock.  Expressing 
the  derived  functions  in  Woolhouse's  formula  in  terms  of  the 
successive  differences  of  f„,  J],  .../,.  we  obtain  the  formula  in 
question, 

vn  -  \ 

/o  +f  1  +/i  +  •  •  •  +/;■  =  m  (f,  -f/i  + . . .  +/.) ~  if, + /;.) 


3.i0        1 

E.rample.  1. — 'I'o  determine       2     ~~ • 

«  =  300     " 

Let   m  =  10,  w  =  10. 

Then  from  Woolhouse's  formvila  we  have 

n  =  zm    '"■   ~    ^300     310     320     330     340     350/ 

2  V3OO    350/      12  V     350-    300V 
=  018512761- 0-02785714 
2432 
=  015724615 

Using  Lubbock's  formula  we  get 
350      1 

-99   ■   1899  (53^-       53/  s  _  99   .  899    t%U-  +  54/M 
-20000      ^    ■'\     "''4'         480000    ^"/s  +  O./s) 

=r0-18512761 -002785714 

1938 

487 

3 

3 

=  015724616 
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To  test  the  accuracy  of  these  answers,  the  reciprocals  of  the  numbers 
300,  301,  302,  ...,  SoO  were  taken  from  Barlow's  tables,  and  summed  by 
means  of  a  comptometer.  The  result  was  found  to  be  0'li'372461G.  The 
closeness  of  the  approximation  is  remarkable. 

176        1 
^a-amp/e  2. —Obtain  the  value  of     2     — • 
n  =  100    "' 

Examph  3. — By  giving  v  increasing  positive  values,  show  that 

T  1 

-L'    (l-t-i  +  ^+  •■■  +  — -log.  «)  =  Or)7721566... 


30.    Formulae  of  Approximate  Quadrature 

Tlie  expan.siojis  whicli  have  been  obtained  in  the  preceding 
articles  furnish  the  value  of  an  integral  to  any  required  degree  of 
accuracy,  provided  a  suHicient  number  of  terms  is  taken.  We  shall 
now  consider  certain  foi"mulae  which  are  in  theii-  nature  only 
approximate,  and  which,  though  frequently  useful,  cannot  be  used 
u  here  great  accuracy  is  required. 

Suppose  that  it  is  required  to  integrate  a  function  J\x)  between 
limits,  which  may,  without  loss  of  generality,  be  taken  to  be  -  1 
and  +  1.     We  shall  endeavour  to  obtain  expressions  of  the  type 

U,f{k)  +  lIJ(h,)+...+IIJ{h„) 

which  represent  the  integral  as  closeh'  as  possible,  where  h^,  h^,  ...  h„ 
are  (n  +  l)  values  of  x  within  the  range  of  integration  and  /<„,  A,,  ...  k,„ 
//,„  //,,  ...  //„  do  not  depend  on  the  function _/"(,i';. 

Let  us  first  choose  //q,  U^,  ...  H,„  so  as  to  make  the  formula 
strictly  accurate,  so  long  as  /(x)  is  a  polynomial  of  degree  less 
than  n  +  \. 

Let  F (x)  =  (x  -  h,)  (x-h^)  ...  {x  -  h„). 

F{x) 
Then  —  is  a  polynomial  of  degree  less  than  7i  +  l,  .so  that  tlie 


fiirmula  must  be  strictly  accurate  when   for  ,/'(.'•;   we  put 
This  gives 

J   -  1    X  -  A,.  r-O  X  -  K 


F(x) 


■h 
=  H,F'{K). 


S-2  INTERPOLATION  AND  NUMERICAL  INTEGRATION     Lch.  iv 

Hence,  whatever  h^^,  //„  ...  /i„  may  be,  in  order  that  the  above  con- 
dition may  be  satisfied,  the  values  of  JIo^  11^,  ...  II„  must  be 
given  by 

1     r+i  Fix)  , 

In  practice  it  is  convenient  to  make  the  intervals  between 
successive  values  of  the  argument  equal  to  each  other.  SupposCj 
then,  that  Ao,  ^i,  ...  A„  are  chosen  so  as  to  divide  the  whole  range 
of  integration  into  ti  equal  parts. 

Then  h,  =    -  I 

A,  =    -1+  — 


/io    =      -    1   + 


H,. 


{h,  -K)  ...  {h,  -  /i,_i ){h,-  /i,+, )  . . .  {h,  -  h„ ) 

X  I      (a:  -  Ao)  . •  •  {x  -  A,._i)  {x  -  h,._^_^)  ...  (x  -  /t„)  di 


2r^    2(r-l)       ^    zl    zA        -  ^  (n  -  r) 
n  n  n        n        ti      "  n 


|>.,)...(...-i<r^')(..,-?(^ti))..,..-,, 

(  -  1)'-^ 


.|:;......(...-^i^>)(,..,-^i^>)..,„.-.). 


Now,  let  t  =  ~^  (x  +  1 ). 


Then 


H,.^ -^ -r       t{t-\)...  {t-r+\){t-r~l)  ...  (t-7i)di 

71 .  r !  (ri-r)  !  Jo  '  '        ^         ' 


Hence,  finally,  if  we  put 

Aq  =  a,  ...  h^  =  a-\-r 


dt. 
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wo  oV)taiii 

J  a  r=U 

wlieie 

J^r--   -r, ,  \    t{t.-\)  ...(t  -r+\){(-r-\)  ...{f.~n) 

r  !  {u  -  r) : Jo 

This  is  known  as  Cotes  Formula  of  Quadrature.  By  means  of 
it  we  may  calculate  the  value  of  the  integral  without  first  of  all 
forming  a  table  of  differences. 

Special  Cases. 

(i.)     Let«=l.     Tlien 

Ca-\  w 

/(x)  dx  =  Hj\a)  +  HJ{a  +  ui) 

-  w     fi 
where  H,,  =  {i,-\)dt  =  },  w 

U  :     1   .   Jo 

Hence 

This  is  known  as  the  Trapezoidal  Rule. 
(ii.)     Letn  =  2.     Then 

r^"  ""/(ic)  dx  =  HJ(a)  +  HJin  4-  w)  +  Iljia  +  2  ?<;) 
whore 

//.=  jjfj,  £('-!)  ('-2) -i'-^c  . 

Hence 

f{x)dx  =  -3-/(«)  +  —^/{n  +  n-)  +  y/(a  +  2  7r). 
This  is  called  Simpson's  *  Hiile  or  the  Parabolic  Rule. 
*  It  was,  however,  first  given  by  James  Gregory. 


84  INTERPOLATION  AND  NUMERICAL  INTEGRATION      [cii.  iv 

Corollary. — If  we  divide  the  whole  range  of  integration,  say 
{a,  b),  into  2  n  equal  parts,  and  apply  Simpson's  Rule  to  each  of 
them,  the  value  of  the  integral  may  be  expressed  as 

/  =  -^  {yn  +  y-2n  +  2  (2/0  +  2/4  +  .  •■  +  2/2„_o)  +  4  (y,  +  y,  +  ...  +  y,„_i)} 
where  y^,  yj,  ...  y._.„  are  the  values  oi  /{x)  when  x  has  the  values 
a,   a  +  — — ,    ...    b.       This   is   Simpson's   Rule  in   its  generalised 

2  71 

form.     The  rule  may  be  stated  thus  : — 

Divide  the  area  into  an  even  number  of  strips  by  equidistant 
ordinates  yo,  yi,  ...  y.,,,'-  t^  ^he  sum  of  the  extreme  orJinates  add 
twice  the  sum  of  the  other  ordinates  with  even  suffixes,  and  four 
times  the  sum  of  the  ordiiiates  with  odd  suffixes,  and  multiply  this 
total  by  one-third  of  the  common  distance  between  the  ordinates* 

(iii)     Let  n  =  3.     Then 

r^"^  V(a^)  f^*'  =  ffof{a)  +  HJ{a  +  iv)  +  U.J  {a  +  2  iv)  +  HJ{a  +  3  to) 

J  a 

where  H,  =  ^^,  ^^{t -\)  {t  -  •2){t -^  dt^%w 

^i=l-f^,|[      t{t~^^){t-'i)dt       =lw 

Hence 

r     '/C^)  '''■'^'  =  f  w;  {/(«)  +  3/(a  +  «;)  +  3/(a  +  2  «>)  +  /(a  +  3  w)}, 

which  is  often  termed  the  Three-Eighths  Rule. 

(iv)  Let  n  =  6.  Then  it  may  be  shown,  as  in  the  above 
special  cases,  that 

f{x)dx  =  ^l^{i\fia)  +  -2lQf(a  +  w)  +  27f{a-i-2w) 

+  272/(«  +  3  m;)  +  27 f  (a  +  4  ?/;)  +  216/(a  +  5  w) 
+  41/(a  +  6t«)}. 

*  The  reader  is  warned  that  the  notation  varies  in  the  different  text- 
books, so  that  in  all  cases  an  examination  of  the  notation  used  should  first 
of  all  be  made. 
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Adding 

TTo  ^'/(«  +  3  «;)  =  1 1^  {/(a)  -  6/(a  +  ?r)  +  \o/{a  +  2  «;) 

-  20/{a  +  3w)  +  ]  b/{a  +  i  w)  -  6/ (a  +  5  w) 

+/{a  +  Qw)'i 
we  obtain 

['  ■  '  7(x)  «J  X  +  T--L_  87  (a  +  3w)  =  ^i^  { -t2/(a)  +  2 1 0/{a  +  w) 

+  i2/{a  +  2w)  +  202/ (a  +  3w) 
+  42/(a  +  iw)  +  210/ (a  +  5  wi 
+  i2/(a  +  6w)} 

=  i\  {/(«)  +  5/(«  +  w')  +y(«  +  2  ^^•) 
+  6/(a  +  3M;)+/(a  +  4M;) 
+  5/(a  +  5  2«)  +/(a  +  6  mj)  } . 

Neglecting    the    term    f  Jo- ^V  («  +  3  ««),    which    will    be    fairly 
small,  we  get  Weddle's  Rule  o/  Quadrature  for  6  intervals. 

J{x)dx  =  f^  {/{a)  +  5/(a  +  w)  +/(a  +  2  w)  +  6/(a  +  3  i«) 

'  +/la  +  Aw)  +  5 /"(a  +  5  t«)  +/(a  +  6  rw)  }. 

31.    Comparison    of    the    Accuracy    of    these    Special 
Formulae. 

We  may  compare  the  accuracy  of  these  formulae  by  evaluating,  by  each 

f7  dx 
of  these  methods,  the  definite  integral   I     -.  ,  the  correct  value  of  which  is 

log  J  =  1-94591. 

(i)     Using  the  Trapezoidal  Rule  six  times,  i.e.,  dividing  the  region  into 
six  equal  portions  and  applying  the  rule  to  each  of  them,  we  get 

n  dx 

—   =  yj-{a)  +/(a  +  w)  +/(a  +  2w)+  ...  +/(«  +  5w)  +  i/(a  +  6  w) 

=  2-02143. 

(ii)     Using  the  Parabolic  Rule  repeated  three  times  we  get 
n  d.v 

=  1-95873. 


INTERrOLATION  AND  NUMERICAL  INTEGRATION     [ch.  iv 

(lii)     Applying  the  Three-Eighths  Rule  twice  we  have 
n  dx 


%  {/{«)  +  3/(a  +  w)  +  3/(«  +  2  lb)  +  2/{a  +  3  w)  +  3/(a  f  4  iv) 

^  +Sf(a  +  ow)+J'{a  +  6w)} 


=  1-96607. 
(iv)     Weddle's  Rule  gives 
dx 

=  1-95286. 


i: 


§30. 


None  of  these  gives  the  result  with  any  high  degree  of  accuracy,  the 
errors  being 

(i.)  U -07552 
(ii.)  0-01282 
(iii.)  0-02016 
(iv.)     0-00695 

n  dx 

Example  1. — Evaluate  I     -7^  by  each  of  the  methods  of  § 

Example  2. — Plot  the  curve  y-  =  2a;  +  25  aud  find  the  area  in  square  inches 
between  the  curve,  the  y-axis,  and  the  double  ordinate  at  a;=  12,  the  unit 
being  one  inch. 

Example  3. — The  ordinates  of  the  boundary  of  the  deck  of  a  ship  are 
8,  30,  36,  40,  42,  42,  40,  38,  34.  28,  and  8  feet  respectively,  and  the  common 
distance  between  them  is  21  feet.     Find  the  area  of  the  deck. 

Example  4- — The  velocity  of  a  train  which  starts  from  rest  is  given  by  the 
following  table,  the  time  being  reckoned  in  minutes  from  the  start  and  the 
speed  in  miles  per  hour. 

t  2  4  6  8       .10         12         14        16         18        20 

mjh        10        18        25        29        32        20         11         5  2     at  rest. 

Estimate  the  total  distance  run  in  20  minutes. 


32.    Gauss'  Method  of  Quadrature. 

The  accuracy  of  the  method  given  in  §  30  may  be  increased  if 
we  no  longer  suppose  the  intervals  between  the  successive  values 
Au,  h^,  ...  A„  of  the  argument  to  be  equal. 

Let  us  now  endeavour  to  choose  H^,  H^,  ...  11^^,  Ji^,  h^,  ...  //,„  in 
such  a  way  that  the  formula 

I    ^  /{x)  dx  =  H,j\h,)  +  Hj\h,)  +  ...  +JI„/{h„) 
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shall  be  accurate  for  all  functions  /(x)  of  degree  equal  to  or  less 
than  (2?i+  1). 
Taking 

/{x)  =  F(x)cl>{x) 

where  (f^ix)  denotes  any  polynoniiul  of  degree  <{n+  1).  and  where 

F  (x)  =  {x-  h,)  (x  -  /</)  (X  -  h,)  . . .  (x  -  h„) 

we  see  that  the  condition  to  be  satisfied  may  be  written 


f      F{x)<l>{x)dx  =  0. 


Now  it   is  well  known  that  if  P„+^  (x)  denotes  the  Legend le 
polynomial  *  of  degi'ee  (n+  1),  we  have 

f       P„^,,{x)cf>{x)dx  =  0 

and,  conversely,  this  condition  is  sufficient  to  determine  the 
Legendre  polynomials.  Hence  we  liave  (save  for  a  nnjltiplicative 
constant) 

F{x)  =  P,.+,(x) 

and  therefore  the  quantities  //(„  /*„  ...  h„  must  be  the  roots  of  the 
equation 

P„+,  {X)  =  0. 

The  following  are  the  values  of  the  constants  /)„,  hi,  ...  h^^,  //„,  //,,  ...  H 
in  the  simplest  cases. 


//«-  1 
H,-  1 


(i) 

When«-1. 

ho- 

1 

/',= 

1 

(ii) 

When  7i  =  2. 

/*o- 

-     N^X 

A,  = 

0 

h.,^ 

.'i 

*  Cf.    Whittaker  and   Watson's  Modem  Analysis  or  Byerly'a  Fourier's 
Series  and  Spherical  Harmonica. 
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(iii)     When?i=3. 

h,  =  -  v/(  V  +  ^5  V30)  ^0  =  i  -  sV  v'SO 

/,,  =  -  V(i^  -  ^5  v'SO)  ^1  =  H gV  v/30 

/,,  =      V (?  -  7^5  V30)  ^2  =  4  +  sV  \'30 

A.,  ^       v'(H  s'a  n/30)  ^3  =  i  -  s's  V30 

(iv)     When  n  =  4. 

322  -  13  ./70 

h^=  -  v/(f-BljV70) 

/lo=  0 

/t3=       xA§-^jv^70) 

A,=       V(HslTv^70) 

f+i    (i< 
Example  1. — Determine  the  value  of  I        oTl 

(i)     Using  3  urdinates 


^^0 

900 

H, 

322+13^/70 

900 

H, 

128 
-225 

H, 

322+13V70 

900 

H, 

322-13V70 

—                   C\f\{\ 

f+1  dt 

J-i3+^ 


1,81,5  1 


9    •  3-  s''i        5)    •    3  9    ■  3+  Vf 

=  0-69312165. 
(ii)     Using  5  urdinates 


f+i    dt      322 -13^/70     1 322+13^70    

J_i3+«~         900  ■  3-  v'(iH6ljv^70)^         900  ■3-v/(f 


h  sim 


128  .  J^       322  +  13\/70  . 1 

+  225  '    3    +  900  '3+  N/(i  -  sli  v/70) 

322-I3\/70  I 

+ 


900  3+^/(a  +  Ji,^/70) 

=  0-693147157. 

The  correct  value  is  0-69314718,  so  that  using  3  ordinates  the  error  is  in  the 
5th  decimal  place  :  using  5  ordinates,  in  the  8th  decimal  place. 


Example  2. — Determine  the  values 

(ii) 


f+l    dx 

f+l     dx 
J_,l  +  .^ 
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MISCELLANEOUS   EXAMPLES. 

1.  The  values  of  a  function /(.r)  for  the  values  1(I5U,  l-QGO,  1070,  1080, 
1-090,  1-100  of  the  argument  are  1-25386,  1-26996,  1-28619,  1-30254,  1-31903, 
1-33565.     U8e  Newton's  formula  of  quadrature  to  show  that  the  value  of 

f  1-100 

f{x)dx 
J  1-0.^0 
is  0  06 47-2. 

2.  The  ordinates  of  a  plane  curve  are  of  the  following  lengths,  35,  4-7, 
5-8,  68,  76,  8-1,  8-0,  7-2  feet,  and  they  are  4  feet  apart.  What  is  the  area 
included  between  the  two  end  ordinates?  and  what  is  the  distance  of  the 
centre  of  gravity  of  that  area  (i)  from  the  extreme  left-hand  ordinate,  and 
(ii)  from  the  base  line? 

3.  Approximate  to  the  values  of  the  following  definite  integrals  : — 

(i)         f'   d.rV(3-x  +  a-3) 


(ii) 


(iii) 


(f-x) 
dx 


/(I- Jtjsin2.v) 


4.  Use  (i)  Lubbock's  formula,  (ii)  Woolhouse's  formula,  to  determine  the 
value  of  an  annuity-certain  for  36  years,  interest  being  at  the  rate  of  3%. 

5.  Verify  that  the  area  of  the  curve  y  =  A-^Bx  +  Cx~  +  D3~>  between  the 
limits  x  =  h  and  x=  -h  \s  equal  to  the  product  of  h  into  the  sum  of  the 
ordinates  at  x  =  h/  ^'3  and  x=  -h/jS. 

In  the  case  of  the  curve  i/=f{x)  =  A  +  Bx  +  Cx"  +  Dir'  -V  Ex^  +  Fx''  verify  in 
like  manner  that  the  area  between  x  =  h  and  x=  -h  is  equal  to 

{ of(h  Jl)  +  »f  (0)  +  5/(  -hjj)}  h/9. 

6.  Corresponding  values  of  .'■  and  y  are  given  in  the  following  table,  the 
unit  being  the  inch.  Find  the  volume  of  the  solid  of  revolution  formed  when 
this  curve  revolves  about  the  a--axis.  Find  also  the  centre  of  gravity  of  this 
uniform  solid. 


0 

4 

8 

1-2 

16 

•20 

0-50 

2-58 

4-41 

5-40 

5-10 

0-50 
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7.  A  body  of  200  lbs.  is  acted  on  b}'  a  variable  force  F.  No  other  force 
acts  on  the  body.  AVhen  the  body  has  passed  through  the  distance  x  feet 
the  force  in  pounds  is  as  follows  :  — 


X      0 

0-1 

0-2 

03 

0-4 

Of) 

0-6 

07 

0-8 

0-9 

1  0 

F    20 

21 

21 

20 

19 

185 

1 8  () 

13  5 

9 

4-5 

0 

Determine  the   work  done   upon   the   body   from  a;  =  0  to  a;  =  0-4.       If   the 
velocity  is  0  when  a;  =  0,  what  is  its  value  when  a;  =  0'4  ? 

8.  The  semi-ordinates  in  feet  of  the  deck  plan  of  a  ship  are  respectively 
3,  16-6,  25-5,  28-6,  298,  30,  29-8,  29-5,  28-5,  24-2,  and  6-8,  the  common 
interval  between  them  being  28  feet.     Find  the  area  of  the  deck. 

9.  The  "half"  ordinates  in  feet  of  the  mid  ship  section  of  a  vessel  are 
12-5,  12-8,  12-9,  13,  13,  128,  12-4,  11-8,  10-4,  6-8,  05  respectively,  and  the 
common  interval  is  2  feet.  Find  the  area  of  the  whole  section  and  the 
position  of  its  centre  of  gravity. 

10.  The  area  of  a  ship's  load-water  plane  is  8000  eq.  ft.  ;  the  body 
below  the  load-water  plane  is  divided  into  six  portions  by  equidistant 
horizontal  sections,  3  feet  apart,  whose  areas  in  sq.  ft.  are  respectively 
7600,  7000,  6000,  4500,  2800,  and  250.  Find  (a)  the  displacement  in  tons, 
and  (h)  the  number  of  tons  which  must  be  taken  out  of  the  ship  to  lighten 
her  4  inches. 
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